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Abstract. We study some Riemannian metrics on the space of smooth regu- 
lar curves in the plane, viewed as the orbit space of maps from S 1 to the plane 
modulo the group of diffeomorphisms of S , acting as reparameterizations. 
In particular we investigate the metric for a constant A > 0: 

Gf(h,k):= [ {l + A Kc (6) 2 ){h{8),k(8))\c'{e)\de 
is 1 

where k c is the curvature of the curve c and h, k are normal vector fields to 
c. The term An 2 is a sort of geometric Tikhonov regularization because, for 
A = 0, the geodesic distance between any 2 distinct curves is 0, while for 
A > the distance is always positive. We give some lower bounds for the 
distance function, derive the geodesic equation and the sectional curvature, 
solve the geodesic equation with simple endpoints numerically, and pose some 
open questions. The space has an interesting split personality: among large 
smooth curves, all its sectional curvatures are > 0, while for curves with high 
curvature or perturbations of high frequency, the curvatures are < 0. 



1. Introduction 

This paper arose from the attempt to find the simplest Riemannian metric on 
the space of 2-dimensional 'shapes'. By a shape we mean a compact simply con- 
nected region in the plane whose boundary is a simple closed curve. By requiring 
that the boundary curve has various degrees of smoothness, we get not just one 
space but a whole hierarchy of spaces. All these spaces will include, however, a 
core, namely the space of all shapes with C°° boundary curves. We expect that 
the most natural shape spaces will arise as the completions of this core space in 
some metric hence we take this core as our basic space. Note that it is the orbit 
space 

B e (S\R 2 ) = Emb(S' 1 ,R 2 )/Diff(S' 1 ) 

of the space of all C°° embeddings of S 1 in the plane, under the action by compo- 
sition from the right by diffeomorphisms of the circle. The space Emb(S' 1 ,R 2 ) is 
a smooth manifold, in fact an open subset of the Frechet space C°°(S 1 , 1R 2 ), and 
it is the total space of a smooth principal bundle with base B e (S 1 ,M 2 ) 
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In fact, most of our results carry over to the bigger orbit space of immersions 
mod diffeomorphisms: 



This action is not quite free (see 2.4 and 2.5), hence this orbit space is an orbifold 
(see 2.5) and not quite a manifold. There is the slightly smaller space Imm/ 
(see 2.1) of immersed curves where diffeomorphisms act freely, the total space of 
a principal fiber bundle with a natural connection admitting parallel transport. 
Existence of horizontal curves, however, holds also in the big space Imm (see 2.5) 
which will be one of the weapons in our hunt for geodesies on Bi. 

The second author was led to study the space B e from its relevance to computer 
vision. To understand an image of the world, one needs to identify the most salient 
objects present in this image. In addition to readily quantifiable properties like 
color and area, objects in the world and their projections depicted by 2D images 
possess a 'shape' which is readily used by human observers to distinguish, for 
example, cats from dogs, BMW's from Hondas, etc. In fact people are not puzzled 
by what it means to say two shapes are similar but rather find this a natural 
question. This suggests that we construct, on some crude level, a mental metric 
which can be used to recognize familiar objects by the similarity of their shapes 
and to cluster categories of related objects like cats. Incidentally, immersions also 
arise in vision when a 3D object partially occludes itself from some viewpoint, 
hence its full 2D contour has visible and invisible parts which, together, form an 
immersed curve in the image plane. 

It is a central problem in computer vision to devise algorithms by which com- 
puters can similarly recognize and cluster shapes. Many types of metrics have been 
proposed for this purpose [7]. For example, there are L 1 -type metrics such are the 
area of the symmetric difference of the interiors of two shapes. And there are 
L°°-type metrics such as the Hausdorff metric: the maximum distance of points 
on either shape from the points on the other or of points outside one shape from 
points outside the other. These metrics will come up below, but the starting point 
of this investigation was whether one could use the manifold structure on the space 
of shapes and define an L 2 -type metric by introducing a Ricmannian structure on 
the space. 

Such questions have also arisen in Teichmuller theory and string theory, where 
the so-called Weil-Peterssen metric on the space of shapes (also called the 'universal 
Teichmuller space) has been much studied. In a second part of this paper, we will 
compare our metric to this remarkable (homogeneous!) metric. 

In this paper, we sought the absolutely simplest Riemannian metric that the 
space Bi supports. The most obvious Diff (5 1 )-invariant weak Ricmannian metric 
on the space of immersions is the 7J°-metric: 



BiCS 1 ,!* 2 ) = Imm(S ,1 ,IR 2 )/Diff(5 1 ). 
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where c : S 1 — > R 2 is an embedding defining a point in B e and h, k are vector fields 
along the image curve, defining two tangent vectors to Imm(5 1 ,R 2 ) at c. This in- 
duces a Diff (S' 1 )-invariant weak Riemannian metric on the space of all immersions 
and on Emb(S' 1 , R 2 ), and for the latter space it induces a weak Riemannian metric 
on the base manifold B e . 

Surprisingly, the Riemannian distance defined as the infimum of the arclcngth of 
paths connecting two points in B e (S 1 ,M 2 ) turns out to be 0, see 3.10! This seems 
to be one of the first examples where this purely infinite dimensional phenomenon 
actually appears. 

Motivated by the proof of this result 3.10 we are led to consider the invariant 
Riemannian metric 3.2.6 for a constant A > 0: 



where k c (0) is the curvature of c at c(9). We will argue that this induces a 
reasonable metric on B e (S , R ), as the infimum of arclengths of paths connecting 
distinct points is always positive. Another reason is that the length function 
t : B e (S 1 ,R 2 ) — > R>o has the following Lipschitz estimate 3.3.2 with respect to 
this Riemannian distance: 



In fact, one can bound the Frechet distance between two curves in terms of this 
metric (see 3.5). The completion of the space of smooth curves in this metric 
contains all curves for whose curvature exists weakly as a finite signed measure 
(e.g. piecewise C 2 curves) and is contained in the space of Lipschitz maps from S 1 
to R 2 modulo a suitable equivalence relation, see 2.11. 

The geodesic equation for the metric G A on Emb(S 1 ,R 2 ) and on B e (S 1 ,M. 2 ) 
can be found in 4.1.1: It is a highly non-linear partial differential equation of order 
4 with degenerate symbol, but which nonetheless seems to have a hypoelliptic lin- 
earization. If A — 0, the equation reduces to a non-linear second order hyperbolic 
PDE, which gives a well defined local geodesic spray. For any A, the sectional 
curvature on B e (S 1 ,R 2 ) has an elegant expression which can be found in 4.6.2 
and 4.6.4. It is non-negative if A = and, for general A, becomes strictly negative 
only if the curve has large curvature or the plane section has high frequency. Of 
course we would have liked to solve the problem of existence and uniqueness of 
geodesies for A > 0. We can, however, translate the minimization of path length 
in our metric into an anisotropic Plateau-like problem: In 3.12 we show that a 
curve projects onto a geodesic in B e (S 1 ,R 2 ) if and only if its graph in [0, 1] x R 2 
is a surface with given boundary at {0} x R 2 and {1} x R 2 which is critical for the 
anisotropic area functional 3.12.3. 

In 5.1 we determine the geodesic running through concentric circles and the 
equation for Jacobi vector fields along this geodesic. The solution of the ordinary 
differential equation 5.1.1 describing this geodesic can be written in terms of elliptic 
functions. This geodesic is no longer globally minimizing when the radius of the 
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circles is large compared to \f~A and has conjugate points when it hits this positive 
curvature zone. In 5.2 we study geodesies connecting arbitrary distant curves, 
hence requiring long translations. The middle part of such geodesies appear to 
be approximated by a uniformly translating 'cigar'-like curve with semi-circular 
ends of radius \f~A connected by straight line segments parallel to the direction of 
translation. These figures were found by numerically minimizing a discrete form 
of the energy functional 3.12.1. 

Finally, in 5.3 and 5.4, we have some further pictures of geodesies. First we 
examine the formation of singularities when a small perturbation is propagated 
forward and A = 0. Then we look at some geodesic triangles in B e whose vertices 
are ellipses with the same eccentricity and center but different orientations. For 
various values of A, we find that these triangles have angle sums greater and less 
than 7r. 

2. The manifold of immersed closed curves 

2.1. Conventions. It is often convenient to use the identification R 2 = C, giving 
us: 

xy = (x,y) +idet(x,y), det(x,y) = (ix,y). 
We shall use the following spaces of C°° (smooth) diffeomorphisms and curves, 
and we give the shorthand and the full name: 

Diff(S' 1 ), the regular Lie group ([6], 38.4) of all diffeomorphisms S 1 — > S 1 with 
its connected components Diff + (5' 1 ) of orientation preserving diffeomor- 
phisms and Diff~(S' 1 ) of orientation reversing diffeomorphisms. 

Diffi(S' 1 ), the subgroup of diffeomorphisms fixing 1 G S 1 . We have diffeomorphi- 
cally Diff(S* 1 ) = Diff^S 1 ) x S 1 = Diff+(5' 1 ) x (S 1 x Z 2 ). 

Emb = Emb(5 1 , R 2 ), the manifold of all smooth embeddings S 1 — > R 2 . Its tangent 
bundle is given by TEmb(S* 1 ,M 2 ) = Emb(S' 1 ,M 2 ) x C 00 ^ 1 , R 2 ). 

Imm = Imm(5 1 ,R 2 ), the manifold of all smooth immersions S 1 — > R 2 . Its tangent 
bundle is given by TImm(5 1 ,R 2 ) = Imm(S' 1 ,M 2 ) x C 00 ^ 1 , R 2 ). 

Imm/ = Immj(S' 1 , R 2 ), the manifold of all smooth free immersions S 1 —> R 2 , 
i.e., those with trivial isotropy group for the right action of Diff (S 1 ) on 
Imm(5 1 ,M 2 ). 

B e = BeiS 1 ^ 2 ) = Emb(S' 1 ,M 2 )/Diff(S' 1 ), the manifold of 1-dimcnsional con- 
nected submanifolds of R 2 , see 2.3. 

Bi = B^S 1 ,®?) = Imm(S' 1 ,R 2 )/Diff(S* 1 ), an infinite dimensional 'orbifold'; its 
points are, roughly speaking, smooth curves with crossings and multi- 
plicities, see 2.5. 

B u = B hf (S\R 2 ) = Imm / (5 1 ,R 2 )/Diff(S' 1 ), a manifold, the base of a principal 
fiber bundle, see 2.4.3. 
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We want to avoid referring to a path in our infinite dimensional spaces like Imm 
or B e as a curve, because it is then a 'curve of curves' and confusion arises when 
you refer to a curve. So we will always talk of paths in the infinite dimensional 
spaces, not curves. Curves will be in E 2 . Moreover, if t [9 c(t,0)) is a 
path, its t-th curve will be denoted by c(t) — c(t, ). By c t we shall denote the 
derivative d t c, and eg = dgc. 



2.2. Length and curvature on Imm(S' 1 ,R 2 ). The volume form on S 1 induced 
by c is given by 

(1) vol : Emb(S' 1 ,E 2 ) -» ^(S 1 ), vol(c) = \c g \d8 
and its derivative is 

(2) dvol(c)(ft) = ( h ? ,C ^ d6. 
We shall also use the normal unit field 

. c B 
n c = i—. 

M 

The length function is given by 

(3) £:lmm(S 1 ,R 2 )->R, £(c) = / \c e \d6 
and its differential is 

(4) dl(c)(h) = [ { -^d6 = [ (ft, ^ - { -^4 L c e )d9 

Js 1 M Js 1 x \ c 0\ M 1 

= — (ft, K(c).icg)d9 = — / (ft, n c )n(c) vol(c) 
Js 1 Js 1 

The curvature mapping is given by 

(5) k : Im m (5M 2 ) - C~(S\R), «(c) = det \ c ^ = <!^£> 

\cg\ 6 \ce\ 6 

and is equivariant so that k(co /) = ±k(c) o / for / e Diff ± (S' 1 ). Its derivative is 
given by 

/ c x , / w,n (ihg,Cgg) (icg.hgg) ,{hg,Cg) 

With some work, this can be shown to equal: 

,„x , / w M yVg) , (Meg) 2 1/1 , (ft, ic e ) , \ 

(7) ck(c)(ft) = - — r^-KeH j — ;— ^ +] — rh — r(— ; — J—jg) ■ 

\cg\ 2 \cg\ \cg\\\cg\ y \cg\ ,6 )e 

To verify this, note that both the left and right hand side are equivariant with re- 
spect to Diff(5 1 ), hence it suffices to check it for constant speed parametrizations, 
i.e. \cg\ is constant and egg = n\cg\icg. By linearity, it is enough to take the 2 cases 
ft = aicg and ft = beg. Substituting these into formulas (6) and (7), the result is 
straightforward. 
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2.3. The principal bundle of embeddings Emb(S' 1 ,R 2 ). We recall some basic 
results whose proof can be found in [6]: 

(A) The set Emb(S' 1 , M?) of all smooth embeddings S 1 — > R 2 is an open subset 
of the Frechet space C°°(S 1 ,M 2 ) of all smooth mappings S 1 — > R 2 with the C°°- 
topology. It is the total space of a smooth principal bundle ir : Emb(S' , R ) — > 
B e (S 1 ,M. 2 ) with structure group DiS(S 1 ), the smooth regular Lie group group of 
all diffeomorphisms of S , whose base B e (S , R ) is the smooth Frechet manifold 
of all submanifolds o/R 2 of type S , i.e., the smooth manifold of all simple closed 
curves in R 2 . ([6], 44.1) 

(B) This principal bundle admits a smooth principal connection described by the 
horizontal bundle whose fiber Af c over c consists of all vector fields h along c such 
that (h,cg) = 0. The parallel transport for this connection exists and is smooth. 
([6], 39.1 and 43.1) 

See 2.4.3 for a sketch of proof of the first part in a slightly more general situation. 
See also 3.2.2 and 3.2.3 for the horizontal bundle Af c . Here we want to sketch the 
use of the second part. Suppose that t (Q c(t,8)) is a path in Emb(S' 1 ,M 2 ). 
Then tt o c is a smooth path in B e (S 1 ,R 2 ). Parallel transport over it with initial 
value c(0, •) is a now a path / in Emb(S' 1 ,R 2 ) which is horizontal, i.e., we have 
(ft, fe) = 0. This argument will play an important role below. In 2.5 below we 
will prove this property for general immersions. 

2.4. Free immersions. The manifold Imm(5' 1 ,R 2 ) of all immersions S 1 — > M 2 
is an open set in the manifold C°°(S 1 ,M 2 ) and thus itself a smooth manifold. An 
immersion c : S* 1 — > R 2 is called free if Diff (S 1 ) acts freely on it, i.e., c o ip = c for 
<p € Diff(5' 1 ) implies ip = Id. We have the following results: 

(1) If <p G Diff (S 1 ) has a fixed point and if co ip = c for some immersion c then 
ip = Id. This is ([2], 1.3). 

(2) If for c € Imm(S' 1 ,R 2 ) there is a point x e c{S l ) with only one preimage then 
c is a free immersion. This is ([2], 1.4). There exist free immersions without such 
points: Consider a figure eight consisting of two touching ovals, and map S 1 to 
this by first transversing the upper oval 3 times and then the lower oval 2 times. 
This is a free immersion. 



(3) The manifold B i j(S 1 ,M 2 ). ([2], 1.5) The set Imm/(/S rl ,R 2 ) of all free im- 
mersions is open in C°° (S 1 ,M. 2 ) and thus a smooth submanifold. The projection 



onto a Hausdorff smooth manifold is a smooth principal fibration with structure 
group Diff(S' 1 ). By ([6], 39.1 and 43.1) this fibration admits a smooth principal 
connection described by the horizontal bundle with fiber N c consisting of all vector 
fields h along c such that (h, eg) = 0. This connection admits a smooth parallel 
transport over each each smooth curve in the base manifold. 



n : Imm / (S ,1 ,R 2 ) 



Imm/(S' 1 , 




= :B U (S\R 2 ) 



Diff(5 1 ) 
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We might view Imm/(S' 1 ,R 2 ) as the nonlinear Stiefel manifold of parametrized 
curves in R 2 and consequently Bij(S 1 ,R 2 ) as the nonlinear Grassmannian of 
unparametrized simple closed curves. 

Sketch of proof. See also [2] for a slightly different proof with more details. For 
celmm/fS 1 ^ 2 ) and s = (si,s 2 ) G V(c) C C°°{S\WL x S 1 ) consider 

tp c (s) : S 1 -> R 2 , p c (*)(0) = c(s 2 (0)) + *i(s2W).n c (s 2 (e)) 

where V(c) is a C°°-opcn neighborhood of (0, Idgi) in C°°(S' 1 ,R x S* 1 ) chosen in 
such a way that: 

• s 2 e Diff(S' 1 ) for each s e V(c). 

• y c (s) is a free immersion for each s € V(c). 

• For (si, s 2 ) € V(c) and a e Diff we have (si, s 2 o a) e V(c). 

Obviously <^ c (si, S2)°<^ = <Pc{si, S2°a) and s 2 is uniquely determined by f c (si, s 2 ) 
since this is a free immersion. Thus the inverse of ip c is a smooth chart for the man- 
ifold Imm/(5 1 , R 2 ). Moreover, we consider the mapping (which will be important 
in section 4 below) 

V> c : C°°(S\ (-£,£)) -» Imm/^.K 2 ), Q(c) := ^(C 00 ^ 1 , (-£,£))) 

Vc(/)W = c(B) + f(6)n c {0) = <p c (fMsi)(0), 

7T o V : (-e,e)) -» BijiS 1 ,^ 2 ), 

where e is small. Then (an open subset of) V(c) splits diffeomorphically into 

C°°(S\ (-£,£)) x Diff S* 1 

and thus its image under <^ c splits into Q(c) x Diff (S 1 ). So the inverse ofnoip c is 
a smooth chart for Bij(S 1 , R 2 ). That the chart changes induced by the mappings 
93 c and -0 C contructed here are smooth is shown by writing them in terms of 
compositions and projections only and applying the setting of [6]. □ 

2.5. Non free immersions. Any immersion is proper since S 1 is compact and 
thus by ([2], 2.1) the orbit space ^(S^R 2 ) = Imm(5 1 , M 2 )/ Diff(S' 1 ) is Hausdorff. 
Moreover, by ([2], 3.1 and 3.2) for any immersion c the isotropy group Diff(5' 1 ) c 
is a finite cyclic group which acts as group of covering transformations for a finite 
covering q c : S 1 — > S 1 such that c factors over q c to a free immersion c : S 1 — > R 2 
with coq c = c. Thus the subgroup Diffi(S' 1 ) of all diffcomorphisms ip fixing 1 G S 1 
acts freely on Imm(S' 1 ,l 2 ). Moreover, for each c € Imm the submanifold Q(c) 
from the proof of 2.4.3 (dropping the freeness assumption) is a slice in a strong 
sense: 

• Q(c) is invariant under the isotropy group Diff(S' 1 ) c . 

• If Q(c) o ip n Q(c) 7^ for ip e Diff (S 1 ) then ^ is already in the isotropy 
group p e Diff(S' 1 ) c . 

• Q(c) o Diff (5 1 ) is an invariant open neigbourhood of the orbit c o Diff (5* 1 ) 
in Imm(5' 1 ,R 2 ) which admits a smooth retraction r onto the orbit. The 
fiber r~ 1 (co p) equals Q(co p). 
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Note that also the action 

Imm(S\R 2 ) x DifF(S' 1 ) -» Imm(5' 1 ,R 2 ) x Imm(5' 1 ,M 2 ), (c,(p) h-> {c,co<p) 

is proper so that all assumptions and conclusions of Palais' slice theorem [8] hold. 
This results show that the orbit space Bi(S 1 ,M 2 ) has only very simple singular- 
ities of the type of a cone C/{e 2Trfc /™ : < k < n} times a Frechet space. We 
may call the space B i (S 1 7 'R 2 ) an infinite dimensional orbifold. The projection 
7T : Imm(5 1 ,M 2 ) -> B^S 1 ,WL 2 ) = Imm(5 1 , M 2 )/ Diff(S* 1 ) is a submersion off the 
singular points and has only mild singularities at the singular strata. The normal 
bundle N c mentioned in 2.3 is well defined and is a smooth vector subbundle of 
the tangent bundle. We do not have a principal bundle and thus no principal 
connections, but we can prove the main consequence, the existence of horizontal 
paths, directly: 

Proposition. For any smooth path c in Imm(5' 1 ,R 2 ) there exists a smooth path 
ip in Diff(S' 1 ) with tp(0, ) = H51 depending smoothly on c such that the path e 
given by e(t,6) = c(t,(p(t,6)) is horizontal: e t -Leg. 

Proof. Let us write e = c o tp for e(t, 9) = c(t, ip(t, 9)), etc. We look for tp as the 
integral curve of a time dependent vector field £(t, 9) on 5 1 , given by tp t = £ <P- 
We want the following expression to vanish: 

(d t (co ip),d g (co <p)) = (c t o ip + (c g o tp) p t , ( Cg o tp) tp e ) 

= ((c t ,c g ) o ip) ipg + {(cg,Cg) o tp) ipg tp t 
= (((c t ,Cg) + (eg, Cg) £) o tp) ipg. 

Using the time dependent vector field ^ = — an d its flow ip achieves this. □ 

2.6. The manifold of immersions with constant speed. Let Imm a (5 1 ,R 2 ) 
be the space of all immersions c : S* 1 — > R 2 which are parametrized by scaled arc 
length, so that \cg\ is constant. 

Proposition. The space Imm a (S' 1 ,R 2 ) is a smooth manifold. There is a diffeo- 
morphism Imm(S' 1 ,R 2 ) = Imm a (S' 1 ,R 2 ) x Diff 1 _ (5' 1 ) which respects the splitting 
Diff (S 1 ) = Difff (S 1 ) x (S 1 x Z2). There is a smooth action of the rotation and re- 
flection group S 1 xZ 2 on Imm a (S' 1 , R 2 ) with orbit space Imm a (5' 1 , R 2 )/(5 1 xZ 2 ) = 
B t (S\R 2 ). 

Proof. For c e Imm(S' 1 ,R 2 ) we put 

, is , , /27ri L \c / (u)\du\ 

^GDiff!^ 1 , <7 c g=exp( / \\ 

a : Imm(5 1 ,R 2 ) -» Imm a (5 1 ,R 2 ), a(c)(9) := c(a~ 1 (9)). 

By the fundamentals of manifolds of mappings [6] the mapping a is smooth from 
Imm(5 1 ,R 2 ) into itself and we have ao a — a. 
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Now we show that Imm a (5 1 ,R 2 ) is a manifold. We use the notation from the 
proof of 2.4.3 with the freeness assumption dropped. For c £ Imm a (5 1 ,R 2 ) we 
use the following mapping as the inverse of a chart: 



The chart changes are smooth: If for (fi,9i) £ C°°(S 1 , (— e, e)) x S 1 we have 
a(V>ci( +0i)(/i( +0i))) = a(^c 2 ( +e 2 )(/2( +#2))) then the initial points agree 
and both curves are equally oriented so that c\{9 + 9\) + f\{6 + 9i)n Cl (9 + 9{) = 
C2(<p(9) + 62) + f2{y{9) + 92)n C2 (ip(9) + 62) for all 9. From this one can express 
{f 2, 92) smoothly in terms of (/1, 

For the latter assertion one has to show that a smooth path through e\ in Q(ci) 
is mapped to a smooth path in Diffi(S' 1 ). This follows from the finite dimensional 
implicit function theorem. The mapping a is now smooth into Imm (5' 1 , K 2 ) 
and the diffeomorphism Imm(5' 1 ,IR 2 ) — ► Imm a (5 1 ,R 2 ) x Diff^S* 1 ) is given by 
c 1 ► (a(c), a c ) with inverse (e, ip) 1— ► e o (p^ 1 . Only the group S 1 x Z 2 of rotations 
and reflections of S 1 then still acts on Imm a (S' 1 ,R 2 ) with orbit space Bi(S 1 ,M. 2 ). 
The rest is clear. □ 

2.7. Tangent space, length, curvature, and Frenet-Serret formulas on 

Imm a (5 1 ,R 2 ). A smooth curve t i-> c ( ,t) £ ImmfS* 1 ,! 2 ) lies in Imm a (S' 1 ,M 2 ) 
if and only if \dgc\ 2 = \ce\ 2 is constant in 9, i.e., dg\cg\ 2 — 2(cg,cgg) = 0. Thus 
h = d t \ c £ T c Imm(5 1 ,R 2 ) = C 00 ^ 1 ^ 2 ) is tangent to Immjs 1 ^ 2 ) at the foot 
point c if and only if (h$,c$g) + (hgg,cg) = (hg 1 cg)g — 0, i.e., (h$,c$) is constant 
in 9. For c £ Imm a (5' 1 , K 2 ) the volume form is constant in 6> since |c#| = £(c)/2n. 
Thus for the curvature we have 



n^KeDxS 1 ^ |J Q(c( +fl))Almm a (5 1 ,R 2 ), 



(/,0)~V>c( +*)(/( +fl))^a(^( +^))) 




and for the derivative of the length function we get 





The derivative of the curvature thus becomes: 
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2.8. Horizontality on Imm a (S' 1 , M 2 ). Let us denote by Imm ai /(S' 1 , R 2 ) the 
splitting submanifold of Imm consisting of all constant speed free immersions. 
From 2.6 and 2.4.3 we conclude that the projection Imm 0)/ (S 1 , R 2 ) -> B f (S 1 ,M 2 ) 
is principal fiber bundle with structure group S 1 x Z 2 , and it is a reduction of 
the principal fibration Imm/ —* Bj. The principal connection described in 2.4.3 
is not compatible with this reduction. But we can easily find some principal 
connections. The one we will use is described by the horizontal bundle with fiber 
N a ,c consisting of all vector fields h along c such that {hg,cg)e = (tangent to 
Imm ) and (h(l),cg(l)) = for 1 <G S 1 (horizontality). This connection admits 
a smooth parallel transport; but we can even do better, beyond the principal 
bundle, in the following proposition whose proof is similar and simpler than that 
of proposition 2.5. 

Proposition. For any smooth path c in Imm (S' 1 , R 2 ) there exists a smooth curve 
(p c in S 1 with <£ c (0) = 1 depending smoothly on c such that the path e given by 
e(t,0) — c(t,ip c (t)6) is horizontal: e t (l)-Leg(l). □ 

2.9. The degree of immersions. Recall that the degree of an immersion c : 
S 1 — > R 2 is the winding number with respect to of the tangent d : S 1 — ► R 2 . 
Since this is invariant under isotopies of immersions, the manifold Imm(5' 1 ,IR 2 ) 
decomposes into the disjoint union of the open submanifolds Imm fc (S\R 2 ) for 
k € Z according to the degree k. We shall also need the space Imm^ 1 , R 2 ) of all 
immersions of degree k with constant speed. 

2.10. Theorem. 

(1) The manifold Imm fc (S' 1 ,R 2 ) of immersed curves of degree k contains the 
subspace Imm*(5' 1 , R 2 ) as smooth strong deformation retract. 

(2) For k =/: the manifold Imm^(5' 1 , M. 2 ) of immersed constant speed curves 
of degree k contains S 1 as a strong smooth deformation retract. 

(3) For k ^0 the manifold B^{S\R 2 ) := Imm*(5 1 , R 2 )/ Diff+(5' 1 ) is con- 
tractible. 

Note that for k ^ Imm' is invariant under the action of the group Diff + (5' 1 ) 
of orientation preserving diffcomorphism only, and that any orientation reversing 
diffcomorphism maps Imm' to Imm - '. 

The nontrivial S 1 in Imm fe appears in 2 ways: (a) by rotating each curve around 
c(0) so that c'(0) rotates. And (b) also by acting S 1 3 (3 i-> (c(6») ^ c(/30)). The 
two corresponding elements a and b in the fundamental group are then related by 
a k = b which explains our failure to describe the topological type of B®. 

Proof. (1) is a consequence of 2.6 since Difff (S 1 ) is contractible. 

The general proof is inspired by the proof of the Whitney- Graustein theorem, 
[9], [4], [3]. We shall view curves here as 27r-periodic plane- valued functions. For 
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any curve c we consider its center of mass 

1 f 2n 

C(c) = Ccnter(c) := — / c{u) \c'(u)\ du £ R 2 
«(c) Jo 

which is invariant under Diff(S' 1 ). We shall also use a(c) = c'(0)/|c'(0)|. 

The case k ^ 0. We first embedd S 1 into Imm(S' 1 , R 2 ) in the following way. For 
a £ S 1 C C = R 2 and fc^Owe but e a {6) = a.e lke /ik, a circle of radius l/\k\ 
transversed fc-times in the direction indicated by the sign of k. Note that we have 
Center(e Q ) = and e' a (0) = a. 

Since the isotopies to be constructed later will destroy the property of having 
constant speed, we shall first construct a smooth deformation retraction A : [0, 1] x 
Imm fe — ► Imm 10 onto the subspace Imm 10 of unit speed degree k ^ curves with 
center 0. 

Let c : R — > R 2 be an arbitrary constant speed immersion of degree k, period 
2ir, and length £(c). Let s c (v) = J Q V \c'(u)\ du be the arc-length function of c and 
put 

A(c, t, u)=(l-t + t^) . (c((l - t)u + t.8- 1 (^u)) t.C(cj) . 

Then A c is an isotopy between c and ci := A(c, 1, ) depending smoothly on c. 
The immersion ci has unit speed, length 27r, and Center(ci) = 0. Moreover, for 
the winding number wo around we have: 

^o(ci|[o,2 W ]) = deg(ci) = deg(c) = k = deg(e Q(c) ) = «Jo(e' a ( c )|[o,27r])- 

Thus Imm fc contains the space Imm 10 of unit speed immersions with center of 
mass and degree k as smooth strong deformation retract. 

For c £ Immf a unit speed immersion with center we now construct an 
isotopy 1 1 — ► H 1 (c, t, ) between c and a suitable curve e a . It will destroy the unit 
speed property, however. For d&rg = - x dy+v dx wc p U ^ : 

y/x 2 +y 2 



if c (u) := [ darg, so that c'(u) = c'(0) e llp " [u) , 

Jc'\[0,u] 
1 /" 27r 

a(c) := — J {ifc(v) - kv)dv, 
ip c (t, u) := (1 - t)ip c {u) + t(ku + a{c)), 

h{c,t,u):= e^'^dv / e^'^dv, 

Jo 27r Jo 

^(c^u) := c'(0)(h(c,t,u) - Center(/i(c, t, fj 

Then H 1 (c, t, u) is smooth in all variables, 27r-periodic in u, with center of mass 
at 0, i? 1 (l,c, u) equals one the e a 's, and i/ 1 (0, c, u) = c(u). But i? 1 (c, t, ) 
is, however, no longer of unit speed in general. And we still have to show that 
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t 1 ► h(c,t, ) (and consequently H 1 ) is an isotopy. 

1 r 2v 

d u h(c, t, u) = e^cC*,") / e^'^dv, 

27r Jo 



(4) 



e 



iV>c(t,«) du 



< 1. 



1 

2^ 

If the last inequality is strict we have d u h(t, u) ^ so that h is an isotopy. If 
we have equality then tp c (t,v) is constant in v which leads to a contradiction as 
follows: If k 7^ then ip c (t, 2ir) — ip c {t, 0) = 2irk so it cannot be constant for any t. 

Let us finally check how this construction depends on the choice of the base 
point c(0). We have: 

<fi c (p+ )(u) = <Pc(l3 + u)-(p c (P), 

a(c(0+ ))=a(c) + k0-<p c (j3), 

ip c{f3+ ) (t,u) = ip c (t,u + (3)-<p c ((3), 

h( C (/3+ ), t, u) = e -'^(/3) (fc( C) ti/3 + „) _ /l ( C; t) /3 )) j 

iJ 1 (c(/3+ ),t,u)= J ff 1 (c,i,/3 + w). 

Let us now deform H 1 back into Immf . For c e Imrrq we consider 

H 2 (c,t,u):=i4(l, J ff 1 (c,t ) ),u), 
H 3 (c,t,u) := H 2 (c 7 t,u + (p H 2 (c) (t)), 

where the iff for a unit speed path / is from proposition 2.8, so that iJ 3 (c) is a hor- 
izontal path of unit speed curves of length 2tt, (i.e., d t H 3 (c, t, 0)_L9 tt |o-ff 3 (c, t, u)). 



The isotopy A reacts in a complicated way to rotations of the parameter, 
but we have A(c(f3+ ), l,u) = A(c, l,^s c (/3) + u). Thus H 3 (c( +(3),t,u) = 
H 3 (c, t,u + (3), so H 3 is equivariant under the rotation group S* 1 C Diff (S 1 ). For 
fc ^ Owe get an equivariant smooth strong deformation retract within ImmJ onto 
the subset {e a : a £ S 1 } C Imuq which is invariant under the rotation group 
S 1 C Diff(S' 1 ). It factors to a smooth contraction on Bf. This proves assertions 
(2 ) and (3 ) for k ^ 0. □ 



2.11 Bigger spaces of 'immersed' curves. We want to introduce a larger 
space containing B^S 1 ,M. 2 ), which is complete in a suitable metric. This will 
serve as an ambient space which will contain the completion of Bi(S 1 ,M. 2 ). Let 
Cont(5' 1 , K 2 ) be the space of all continuous functions c : S 1 —> M 2 . Instead of a 
group operation and its associated orbit space, we introduce an equivalence relation 
on Cont(5 1 ,K 2 ). Define a subset R C S 1 x S 1 to be a monotone correspondence 
if it is the image of a map 

x — > (h(x) mod 2-7T, k{x) mod 27r), where 

h, k : M — > R are monotone non-decreasing continuous functions such that 
h(x + 2tt) = h(x) + 2tt, k(x + 2tt) = k(x) + 2tt. 



RIEMANNIAN GEOMETRIES ON SPACES OF PLANE CURVES 



13 



In words, this is an orientation preserving homcomorphism from S 1 to S 1 which 
is allowed to have intervals where one or the other variable remains constant 
while the other continues to increase. (These correspondences arise naturally in 
computer vision in comparing the images seen by the right and left eyes, see [1].) 
Then we define the equivalence relation on Cont(5' 1 , K 2 ) by c <~ d if and only if 
there is a monotone correspondence R such that for all d,tp e R,c(6) — d(ip). 
It is easily seen that any non-constant c e Cont(5' 1 ,E 2 ) is equivalent to an ci 
which is not constant on any intervals in S 1 and that for such ci's and di's, the 
equivalence relation amounts to c\ o h = d\ for some homeomorphism h of S 1 . Let 
Bf ont (S 1 ,R 2 ) be the quotient space by this equivalence relation. We call these 
Frechet curves. 

The quotient metric on Bf ont (S 1 , R 2 ) is called the Frechet metric, a variant of 
the Hausdorff metric mentioned in the Introduction, both being L°° type metrics. 
Namely, define 



It is straightforward to check that this makes Bf ont (S 1 , R 2 ) into a complete metric 
space. 

Another very natural space is the subset Bf p (S 1 ,R 2 ) C Bf ont (S 1 ,M 2 ) given by 
the non-constant Lipschitz maps c : S 1 — > R 2 . The great virtue of Lipschitz maps 
is that their images are rectifiablc curves and thus each of them is equivalent to 
a map d in which 9 is proportional to arclength, as in the previous section. More 
precisely, if c is Lipschitz, then eg exists almost everywhere and is bounded and 
we can reparametrize by: 



obtaining an equivalent d for which \dg\ = L/2n. This d will be unique up to 
rotations, i.e. the action of S 1 in the previous section. 

This subspace of rectifiable Frechet curves is the subject of a nice compactness 
theorem due to Hilbert, namely that the set of all such curves in a closed bounded 
subset of R 2 and whose length is bounded is compact in the Frechet metric. This 
can be seen as follows: we can lift all such curves to specific Lipschitz maps c whose 
Lipschitz constants are bounded. This set is an equicontinuous set of functions by 
the bound on the Lipschitz constant. By the Ascoli-Arzela theorem the topology 
of pointwise convergence equals then the topology of uniform convergence on S 1 . 
So this set is a closed subset in a product of S 1 copies of a large ball in R 2 ; this 
product is compact. The Frechet metric is coarser than the uniform metric, so our 
set is also compact. 




inf 

homcomor ph .h:S 1 ^S 1 



\\coh-d\\ 
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3. Metrics on spaces of curves 



3.1 Need for invariance under reparametrization. The pointwise metric on 
the space of immersions Imm(S' 1 ,R 2 ) is given by 



This Riemannian metric is not invariant under reparameterizations of the variable 
9 and thus does not induce a sensible metric on the quotient space Bi(S 1 ,M. 2 ). 
Indeed, it induces the zero metric since for any two curves Cq,C\ G Bi(S 1 ,M 2 ) 
the inftmum of the arc lengths of curves in Imm(S' 1 ,R 2 ) which connect embeddinqs 
Co,Ci € Imm(5 1 ,R 2 ) with 7r(c,) = d turns out to be zero. To see this, take any 
Co in the Diff (S' 1 )-orbit over Co- Take the following variation c(6,t) of c$: for 9 
outside a small neighborhood U of length e of 1 in S 1 , c(0,t) = cq(9). If 9 € U, 
then the variation for t £ [0, 1/2] moves the small part of Co so that c(6, 1/2) for 9 
in U takes off Co, goes to Ci, traverses nearly all of Ci, and returns to Co- Now 
in the orbit through c(-, 1/2), reparameterize in such a way that the new curve is 
diligently traversing C\ for 9 £ U, and for 9 e U it travels back to Co, runs along 
Co, and comes back to C\. This reparametrized curve is then varied for t € [1/2, 1] 
in such a way, that the part for 9 <E U is moved towards C2. It is clear that the 
length of both variations is bounded by a constant (depending on the distance 
between Co and C\ and the lengths of both Co and C\) times e. 

3.2. The simplest Riemannian metric on Bi. Let h,k <G C°° (S 1 ,R 2 ) be two 
tangent vectors with foot point c <G Imm(5 1 ,lR 2 ). The induced volume form is 
vol(c) = {dec, dgc) 1 / 2 d9 = \cg\d9. We consider first the simple H° weak Riemann- 
ian metric on Imm(5 1 ,R 2 ): 



which is invariant under Diff(S' 1 ). This makes the map ir : Imm(5 1 ,R 2 ) — > 
Bi(S , K ) into a Riemannian submersion (off the singularities of Bi(S , R )) 
which is very convenient. We call this the H° '-metric. 

Now we can determine the bundle TV — > Imm(5 1 ,R 2 ) of tangent vectors which 
are normal to the Diff (5' 1 )-orbits. The tangent vectors to the orbits are T c (c o 
Diff(5' 1 )) = {g.cg : g e C°°(5 1 ,R)}. Inserting this for k into the expression (1) of 
the metric we sec that 




(1) 




(2) 



Af c = {h e C°°(5 1 ,R 2 ) : (h, c e ) = 0} 

= {aicg e C°°(5 1 ,R 2 ) : a E C°°(5 1 ,R)} 
= {bn c G C°°(S\R 2 ) : b G C°°(S' 1 , M)}, 



where n c is the normal unit field along c. 
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A tangent vector ft e T c Imm(5 1 , M 2 ) = C°° (S 1 , M 2 ) has an orthonormal de- 
composition 

h = h T +h ± e T c (coDiff+(5 1 ))eAA c where 

(3) h T = i^l Cee r c ( C oDiff + (^ 1 )), 

± (ft, ice) . ^ . r 
ft = -] — n^-zce G A/c, 

into smooth tangential and normal components. 

Since the Riemannian metric G on Imm(S' 1 ,R 2 ) is invariant under the action 
of Diff(S' 1 ) it induces a metric on the quotient Bi(S 1 ,M. 2 ) as follows. For any 
Cq,Ci G Bi, consider all liftings co,ci G Imm such that 7r(co) = (7o,7r(ci) = C\ 
and all smooth curves t ^ {9 ^ c(t,6)) in Imm(S' 1 ,R 2 ) with c(0, •) = c and 
c(l, •) = Ci. Since the metric G is invariant under the action of Diff(S' 1 ) the 
arc-length of the curve t 1— > n(c(t, •)) in _B i (5 ,1 ,K 2 ) is given by 

L^ OT (c) := L G (n(c(t, •))) - ^ ^ G^ c) {T c iT.c tl T c iT.c t ) dt = J JcJ^f)dt 

f 1 ( [ I (ct,ice) . (c t ,ic g ) . \. ,,^2,, 




The metric on B^S 4 , M 2 ) is defined by taking the infimum of this over all paths c 
(and all lifts c , ci): 

dist^(Ci,C 2 ) = ML^ OI (c). 

Unfortunately we will sec below that this metric is too weak: the distance that it 
defines turns out to be identically zero! For this reason, we will mostly study in 
this paper a family of stronger metrics. These are obtained by the most minimal 
change in G. We want to preserve two simple properties of the metric: that it 
is local and that it has no derivatives in it. The standard way to strengthen the 
metric is go from an H° metric to an H 1 metric. But when we work out the 
natural H 1 metric, picking out those terms which are local and do not involve 
derivatives leads us to our chosen metric. 

We consider next the H 1 weak Riemannian metric on Imm(5 1 ,l 2 ): 

(5) Gl(h,k):= f ((h(9),k(9))+A { -^^)\c e \de. 

Js 1 \ c o\ 

which is invariant under Diff(S' 1 ). Thus ir : Imm(5 1 ,R 2 ) -» B^S 1 ^ 2 ) is again 
a Riemannian submersion off the singularities of Bi(S 1 ,M 2 ). We call this the 
H 1 -metric on Bi. 
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To understand this metric better, we assume h = k = + &]f^y- Moreover, 

h_ 



for any function /(#), we write f s = A for the derivative with respect to arc 



length. Then: 



h s = P^-r = (aic s + bc s ) s = (a s + nb)ic s + (b s - na)c s . 



Therefore: 

Gl(h,h)= { (a 2 + b 2 + A(a s + nb) 2 + A(b s - na) 2 )ds 
Js 1 

= I (a 2 (l + AK 2 )+Aa 2 s )+2AK(a s b-b s a) + (b 2 (l + AK 2 ) + Ab 2 s )ds 
Js 1 

Letting T x and T 2 be the differential operators T x = I + An 2 - A(-^) 2 , T 2 = 
A(k s + 2k^), then integrating by parts on 5 1 , we get: 

Gl(h,h) = ( (T 1 (a).a + 2T 2 (a).b + T 1 (b).b)ds. 
Js 1 

Note that T\ is a positive definite self-adjoint operator on functions on c, hence it 
has an inverse given by a Green's function which we write T-j -1 . Completing the 
square and using that T± is self-adjoint, wc simplify the metric to: 

Gl(h, ^) = ^(T 1 (a). a -T 1 - 1 (r 2 (a)).T 2 (a)+r 1 (6+T 1 - 1 (r 2 ( a ))).(fe+T 1 - 1 (T 2 ( a ))))rf S . 

If we fix a and minimize this in b, we get the bundle Af 1 — ► Imm(S' 1 ,]R 2 ) of 
tangent vectors which are G 1 -normal to the Diff (5' 1 )-orbits. In other words: 

Nl ={heC°°(S 1 ,m 2 ) :h = aic s +bc s ,b= -T^{T 2 {a))} 

and on horizontal vectors of this type: 

G\{h,h)= f((l + Ak 2 )ci 2 + Aa 2 s )ds- / Tf 1 {T 2 {a)).T 2 {a)ds. 

J c J c 

If we drop terms involving a s , say because we assume \a s \ is small, then what 
remains is just the integral of (1 + An 2 )a 2 plus the integral of T^ 1 (n s a)n s a. The 
second is a non-local regular integral operator, so dropping this we are left with 
the main metric of this paper: 

Gc(h,h)= 1(1 + An 2 )a 2 ds,h = aic s 

J c 

which we call the H® -metric with curvature weight A. For further reference, on 
Imm(S' 1 , R 2 ), for a constant A > 0, it is given by 

(6) G£{h,k):= f (l + A Kc (d) 2 )(h(9),k(6))\cJ(9)\de 

Js 1 

which is again invariant under Diff(S' 1 ). Thus ir : Imm(S' 1 ,IR 2 ) — > B^S^R 2 ) is 
again a Riemannian submersion off the singularities. Note that for this metric (6), 
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the bundle J\f C rimm(5' 1 , K 2 ) is the same as for A = 0, as described in (2). The 
arc-length of a curve t 1— > n(c(t, •)) in Bj(5 1 ,K 2 ) is given by the analogon of (4) 

L^J(c) := L GA (n(c(t,-))) = J ^G^ c) (T c n.c t ,T c n.c t ) dt = £ yj dt 



(7) 



(1 + A K 2 c )(c t7 n c ) 2 \c e \dey dt 



The metric on .B^S" 1 , M 2 ) is defined by taking the infimum of this over all paths c 
(and all lifts cq,c\): 

dist^(Ci,C 2 ) = mf L^I(c). 

Note that if a path 7r(c) in B i (S 1 ,R 2 ) is given, then one can choose its lift to a 
path c in Imm(S' 1 ,R 2 ) to have various good properties. Firstly, we can choose the 
lift c(0, ) of the inital curve to have a parametrization of constant speed, i.e. if 
its length is I, then |ce|(0, 0) = £/2ir for all 9 e S 1 . Secondly, we can make the 
tangent vector to c everywhere horizontal, i.e. < c t ,cg >= 0, by 2.5. Thirdly, we 
can reparametrize the coordinate t on the path of length L so that the path is 
traversed at constant speed, i.e. 



(1 + Ak z c ) < c t , ice > dQ/\c e \ = L , for all < t < 1. 



3.3. A Lipschitz bound for arc length in G A . We apply the Cauchy-Schwarz 
inequality to the derivative 2.2.4 of the length function along a path t i— » c(t, ): 



d t £(c) = dt{c){ct) = n(c)(ct,n c )\ce\d0 < / K(c)(c u n c )\c g 

Js 1 Js 1 

< (J l 2 \cg\d6) 5 (J K(c) 2 (c tl n c ) 2 \ce\d9) 5 

< *(c)* ^= (J (1 + An{c) 2 ){c tl n c ) 2 \c e \de) k 



d6 



Thus 



2^1{c) 2y/JyJ S 
and by using (3.2.7) we get 



(^(l + A^c) 2 )^,^) 2 !^)- 



v^i7-^m= r d t {vw))dt 

Jo 



(i) 



U=A L ^ {C) - 
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If we take the infimum over all paths connecting Co with the Diff (S' 1 )-orbit through 
ci we get: 

Lipschitz continuity of v^:B 4 (S ,1 ,M 2 )^R> . For C andd in B^S 1 ,M 2 ) = 
Imm(5 1 ,R 2 )/Diff(5 1 ) we have for A > 0: 

(2) V^Ci)- V^iCo) < ^dist^(Ci,C 2 ). 



3.4. Bounding the area swept by a path in Bi. Secondly, we want to bound 
the area swept out by a path starting from Co to reach any curve C\ nearby 
in our metric. First we use the Cauchy-Schwarz inequality in the Hilbert space 
L 2 {S\\c g (t,9)\d8) to get 

/ i.\ct(t,e)\\c e (t,e)\de = (i,\ Ct \) L 2 < 

Js 1 

< plMMi^ = (J \cg(t,e)\My (j \c t (t,9)\ 2 \c g (t,e)\dey. 

Now we assume that the variation c(t, 6) is horizontal, so that (c t ,cg) = 0. Then 
Lg a { c ) = L^a{c). We use this inequality and then the intermediate value theorem 
of integral calculus to obtain 



L^i{c) = L GA (c) = f y/G£(ct,Ct)dt 

=fS L ° 



(1 + A K (c) 2 )\c t (t,9)\ 2 \cg(t,9)\de] 2 dt 



>-i\; 

Jo vs 



> 



(y \c t (t,8)\ 2 \c e (t,9)\d9ydt 
£ (J \cg(t, 8)\d9y J J \ct(t, 9)\\c e {t, 9)\ d8 dt 
\c e (t ,9)\d9y* I [\ct(t,6)\\cg(t,6)\d9dt 



s 1 ' Jo Js 1 

for some intermediate value < to < 1, 

1 / \detdc(t,8)\d9dt. 



Area swept out bound. If c is any path from Co to C\, then 

( area of the region swept\ < ^ > ~ „ _ jhor,^ 
v ' y out by the variation c J ~ t v v v ' '' G A \ >■ 



This result enables us to compare the double cover Bf r (S 1 ,M. 2 ) of our metric 
space Bi(S-M. 2 ) consisting of oriented unparametrized curves to the fundamental 
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Figure 1. Two distinct immersions of S 1 in the plane whose 
underlying currents are equal. One curve is solid, the other 
dashed. 



space of geometric measure theory. Note that there is a map hi from B° r to the 
space of 1-currents T>[ given by: 

< /ii(c mod Diff+(5 1 )),w >= / c*uj, celmmfS 1 ,! 2 ). 

Js 1 

The image h\(C) is, in fact, closed. For any C, define the integer- valued measur- 
able function wc on R 2 by: 

wc{{x,yj) — winding number of C around (x,y). 

Then it is easy to see that, as currents, h\{C) = d(wcdxdy), hence dh\(C) = 0. 

Although hi is obviously injective on the space B e , it is not injcctivc on Bi as 
illustrated in Figure 1 below. The image of this mapping lies in the basic subset 
Zi iC C T>[ of closed integral currents, namely those which are both closed and 
countable sums of currents defined by Lipschitz mappings Cj : [0, 1] — > R 2 of finite 
total length. Integral currents carry what is called the flat metric, which, for closed 
1-currents, reduces (by the isoperimetric inequality) to the area distance 

(2) dt{C l ,C 2 ) = ff \w Cl -w C2 \dxdy. 

J Jr 2 

To connect this with our 'area swept out bound', note that if we have any path 
c in Imm(5 1 ,R 2 ) joining C\ and C 2 , this path defines a 2-current w(c) such that 
dw(c) = /ii(Ci) - /ii(C 2 ) and 

/ \w(c)\ dxdy < / / \detc\d9dt 
Jk. 2 Jo Js 1 

which is what we are calling the area swept out. But d{wc 1 — wc 2 ) — h\{C\) — 
h\{C2) too, so w(c) — wc 1 — wc 2 - Thus 

(3) d b (C\, C'2) < min [area swept out by c] 

all paths c joining C1.C2 

Finally, we recall the fundamental compactness result of geometric measure 
theory in this simple case: the space of integral 1-currcnts of bounded length is 
compact in the flat metric. This implies that our 'area swept out bound' above 
has the Corollary: 

Corollary. 
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(4) If{C n } is any Cauchy sequence in Bi for the metric dist^A, then {hi(C n )} 
is a Cauchy sequence in Zi iC on which length is bounded. 

(5) Hence hi extends to a continuous map from the completion Bi of Bi in 
the metric G A to I\, c - 



3.5. Bounding how far curves move in small paths in Bi. We want to 
bound the maximum distance a curve Co can move on any path whose length is 
small in G A metric. Fix the initial curve Co and let t be its length. The result is: 

Maximum distance bound. Let e < min{2\/A£, £ 3 ^ 2 }/8 and consider rj = 
4(f 3 / 4 A~ 1 / 4 +^ 1//4 ) v / e. Then for any path c starting at Co whose length is e, the 
final curve lies in the tubular neighborhood of Co of width r\. More precisely, if we 
choose the path c(t, 9) to be horizontal, then maxg |c(0, 9) — c(l, 9)\ < rj. 

Proof. For all of this proof, we assume the path in Bi has been lifted to a 
horizontal path c e Imm(5' 1 ,R 2 ) with \ce\(0, 0) = £/2tt, so that (c t ,c d ) = 0, 
and also J s i(l + AK;?)|c t | 2 |cg| d9 = e 2 . The first step in the proof is to refine 
the Lipschitz bound on the length of a curve to a local estimate. Note that by 
horizontality 




Now we make the key definition: 



MM) = min M(ti,0). 

Note that the t-derivative of \cg\ is either or equal to that of \c$\ and is < 0. 
Thus we have: 




To make use of this inequality, let E = {9 : \c e \(l,6) < (l - (Ae)- 1 / 4 ^/e)£/2Tr}. 
Our assumption on e gives (Ai) -1 ^ 4 ^ < 1/2, hence on S 1 \E we have \cg\ > i/in. 
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On E we have also (M) 1/2 < (1 - (A£)- 1 / 4 ^e/2) v / £/2n. Combining this with 
the previous inequality, we get (where fi(E)is the measure of E): 

We now use the lower bound on |cg| on S 1 — E to control c(l, 9) — c(0, 9): 
f |c(l,0)-c(O,0)|d0< f f \ct\d9dt 

JS 1 -E JO JS 1 -E 

< >/2tt- / if \c t \ 2 d9\ 1/2 dt 

«/ J — E 



V ^7T / / / ,9, , ,/,V' , ZV^TT 



Again, introduce a small exceptional set F = {6 | ^ _E and |c(l,0) — c(0,#)| > 
l 1 / 4 ^}. By the inequality above, we get: 

fj,(F) ■ £ 1/4 y/s < ' hence < < 

The last inequality follows from the second assumption on e. Knowing fJ,(E) and 
/z(-F) gives us the lengths |c(0, and |c(0, F)| in R 2 . But we need the lengths 
\c(l, E) \ and |c(l, F)\ too. We get these using the fact that the whole length of C\ 
can't be too large, by 3.3: 

y/\CT\ <Vf. + -4=, hence 
2\/T 

IT £ 3 / 4 
|Ci| <£ + 2e\ - <£+ v^' 



On we have \c e \ > (1 - (AfT 1 / 4 ^)^-, thus we get 

|c(l,£UF)| = \Cr\-\c{\,S x \{E\lF))\ 

£ 3/4 ( , Vi \ f 



< 



<Vi'(3^ + ^ 1/4 

Finally, we can get from c(0, 6*) to c(l,#) by going via c(0,9') and c(l,#') where 
9' g S 1 \(£UF)^. Thus 

max|c(O,0)-c(l,0)| < |c(0, £ U F)| + £ x l 4 \fe + \c(E UF,1)| 



< 4(^3/4^-1/4 +i ,l/4 ) ^ /i n 



Combining this bound with the Lipschitz continuity of the square root of arc 
length, we get: 
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3.6. Corollary. For any A > 0, the map from Bi(S ,R ) in the dist^ metric to 
the space Bf ont (S 1 , M 2 ) in the Frechet metric is continuous, and, in fact, uniformly 
continuous on every subset where the length I is bounded. In particular, dist G A 
is a separating metric on B i (S 1 ,M. 2 ). Moreover, the completion B^S 1 ,R 2 ) of 
B^S 1 ,^. 2 ) in this metric can be identified with a subset of £?f p (S 1 , R 2 ) . 

If we iterate this bound, then we get the following: 

3.7. Corollary. Consider all paths in Bi joining curves Co and C\. Let L 

be the length of such a path in the dist G A metric and let £ m i n ,£max be the mini- 
mum and maximum of the arc lengths of the curves in this path. Then there are 
parametrizations Cq,C\ of Co and C\ such that: 



To prove this, you need only break up the path into a minimum number of pieces 
for which the maximum distance bound 3.5 holds and add together the estimates 
for each piece. We will only sketch this proof which is straightforward. The 
constant 50 is just what comes out without attempting to optimize the bound. 
The second option for bound, 50^/I^^LF* is just a rephrasing of the bound 
already in the theorem for short paths. If the path is too long to satisfy the 
condition of the theorem, we break the path at intermediate curves d of length 
£i such that each begins a subpath with length Si = min {y/Afi, £ 3 / 2 )/8 and which 
don't overlap for more than 2:1. Thus J2i £ i — Then apply the maximum 
distance bound 3.5 to each piece, letting r/i be the bound on how far points move 
in this subpath or any parts thereof and verify: 



from which we get what we need by summing over i. 

3.8. A final Corollary shows that if we parametrize any path appropriately, we 
get explicit equicontinuous continuity bounds on the parametrization depending 
only on L, £ max and £ m - m . This is a step towards establishing the existence of weak 
geodesies. The idea is this: instead of the horizontal parametrization (c t , eg) = 0, 
we parametrize each curve at constant speed \cg\ = £{t)/2ir where £(t) is the length 
of the t th curve and ask only that (c t , cg)(0, t) = for some base point € [0, 2w], 
see 2.8. Then we get: 

Corollary. // a path c(t, 6), < t < 1 satisfies 



max | c (9) — ci(9)\ < 50max(LF*, y%r^LF*), where 




T]i < 2V2£i < l§V2e t F*, 



\c e (e,t)\=£(t)/27T 
(c t ,c e )(0,t) = 



for all 9, t 
for all t and 
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ther 



MM c(t 2 ,e 2 )\ < ^\e, e 2 \ + 7(e 3 JL/A 1/4 + tULWHt 1 -t 2 ) 



whenever \h - t 2 \ < min(2 v / Al~, £ 3 ^ n )/(8L). 

Proof. We need to compare the constant speed parametrization here with the 
horizontal parametrization - call it c* - used in the maximum distance bound 3.5. 
Under the horizontal parametrization, let the point on Ct 1 correspond to 

{t 2 ,6D on C t2 , i.e. c(t 2 ,^) = c*(t 2 ,0i)- Let C = (iUL/A 1 / 4 + t]lL). Then we 
know from 3.5 that 

|c(ti,0i)-c(t 2 ,0I)| <4CWL(h-t 2 ). 

To compare 9\ and 9{, we use the properties of the set E in the proof of 3.5 to 
estimate: 



(01 - 



2n = / \4(t2,^-^ 

0i£ 2 



> 



2ir 



> - ^(j^y/l 2 - ~ l £ i " ^1 and similarly 
((2ir-0* 1 )-(2ir-9 1 ))£ 2 _ [ 2v . ... 



/ \c* e (t 2 ,y)\dip- 

Jd 1 



2tt J 0l 1 oy r/l r 2tt 

Combining these and using the Lipschitz property of length, we get: 

m ~ ^ l - < 2Cy/L(h -t 2 ) + 2\y/h- y/l 2 \ y/I^ 

< 2C^L(h - t 2 ) + vt^^ < ^Cy/Hh-h) 

VA 2 

Thus, finally: 

|c(ti,0i)-c(i 2 ,0 2 )| < |c(ti,«i) — c(* 2 ,«D| + 

+ | c (t 2 , et) - c(t 2 , e^i + |c(t 2 , ft) - c(t 2 , fl 2 )| 

< 4CVi(ii - t 2 ) + ^CVi(*i - t 2 ) + - 2 |. □ 



3.9. One might also ask whether the maximum distance bound 3.5 can be strength- 
ened to assert that the 1-jets of such curves C must be close to the 1-jets of Co- 
The answer is NO, as is easily seen from looking a small wavelet-type perturbations 
of Co- Specifically, calculate the length of the path: c(t, 9) — cq(9) + t ■ af(9/a) ■ 
*( c o)e(#)jO < t < 1 where f(x) is an arbitrary C 2 function with compact support 
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and a is very small. We claim the length of this path is 0(-Ja), while the 1-jet at 
the point 9 — of the final curve of the path approaches (1 + i/'0)(co)g(0). 

We sketch the proof, which is straightforward. Let C a ,t be the curves on this 
path. Then sup|c t | = O(a), sup\nc at \ = 0(l/a), A < \cg\ < B for suitable 
A,B > and ^(support(c t )) = 0(a). Then the integral J gl (1 + An 2 c ){c tl ice) 2 -^ 
breaks up into 2 pieces, the first being 0(a 2 ), the second being O(l) and the 
integral vanishing outside an interval of length 0(a). Thus the total distance is 
O(v^). 

3.10. The #°-distance on J B 4 (S' 1 ,R 2 ) vanishes. Let c ,ci e lmm(S' 1 ,M 2 ) 
be two immersions, and suppose that t i— > (9 i— > c(t,8)) is a smooth curve in 
Imm(5 1 ,R 2 ) with c(0, •) = c and c(l,-) = a. 

The arc-length for the i?°-metric of the curve t i— ► n(c(t, ■)) in Bj(S' 1 ,K 2 ) is 
given by 3.2.7 as 




Theorem. For Co,Ci € Imm(S ,1 ,R 2 ) there exists always a path t i— » c(t, ■) urci/i 
c(0, •) = Co tmd 7r(c(l, •)) = 7r(ci) smc/i iftai Lg° r (c) is arbitrarily small. 

Heuristically, the reason for this is that if the curve is made to zig-zag wildly, 
say with teeth at an angle a, then the length of the curve goes up by a factor 
1/ cos(a) but the normal component of the motion of the curve goes down by the 
factor cos(a) - and this normal component is squared, hence it dominates. 

Proof. Take a path c(t,9) in Imm(5' 1 ,R 2 ) from Co to c\ and make it horizontal 
using 2.5 so that that (c t , eg) = 0; this forces a reparametrization on c\. 

Now let us view c as a smooth mapping c : [0, f] x [0, f] — * M 2 . We shall use the 
piecewise linear reparameterization ((p(t,0),9) of the square shown above, which 
for < t < 1/2 deforms the straight line into a zig-zag of height 1 and period n/2 
connecting the two end-curves, and then removes the teeth for 1/2 < t < 1. In 
detail: Let c(t,6) = c(<p(t,6),6) where 

'2t(2n9-2k) for < t < 1/2, §f < 9 < ^§±± 

2t(2k + 2-2n9) for < t < 1/2, 2|±i < ff < 2fc±2 



tp(t,0) = { 



2t - 1 + 2(1 - t)(2n9 -2k) for 1/2 < t < 1, || < < ^f±± 

2t - 1 + 2(1 - i)(2fc + 2 - 2n6») for 1/2 < t < 1, ^ < < 2|±2. 
Then we get eg — (pg.c t + eg and c t = <pt-c t where 

' +4nt ( 4n9 - 4fc 

-Ant _ I 4fc + 4 - 4n# 

" ~ ^ +4n(l - t) ' V * ~ | 2 - 4n6> + 4A; 

-4n(l-i) [-(2-4n6» + 4fc) 
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Figure 2. The reparametrization of a path of curves used to 
make its length arbitrarily small. 



Also, (c t ,c e ) = implies (c tl ic g ) = ip t .\c t \.\c e \ and \c g \ = \c e \y/l + ip 2 B {\c t \ / \c g \) 2 . 
Thus 



1 - /•! 



L hor (P 



<^> 2 ^) 



■dt 



1 " 1 rf\ct\ 2 \c e \ , n \\ 



d6_y 

\ce\) Jo yJo ^/l + ^(M )2 



d6 ) dt 



2k 

\k=0 " 2n 



y/ 1 + (4»*) 2 (fe» 



2fc+2 



— (4fc + 4-4nfl) 2 |c t (^,0)| 2 M<p, 



2fc+l 
2n 



>/ 1 + ( 4n *) 9 (fefeffl) a " 

, /„_i 2fc+l 

+ / ( E( /~ (2 ~ 4ng+4fc) |ct( ^ g)| !Mg^W 

"'2 \fe=0 



'-de 



^l + (4n) 2 (l-t) 2 (^) 2 



+ 



2k+2 



— (2 - 4n9 + 4k) 2 \c t (<p,e)\ 2 \c e (<p,0)\ 



2k+l 
2n 



^l + (4n) 2 (l-t) 2 (M^ )2 



d0) dt 



The function \ce(<p, 6)\ is uniformly bounded above and away from 0, and \ct(<p, 6)\ 
is uniformly bounded. Thus we may estimate 

„ i 2fc+l 

g (4nfl~4£Q 2 |c t (^)| 2 M^)| ^ 



f 2fc 

fc=0 " 2n 



v /l + ( 4 ^ 2 (feteST^ 
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0(l) ^i 



. , /M , v 4n^|c t Mt,f + g)|' jfl 



— o + (4nt)2| Ct (^(t, § + 0), § + (9)|2 

We estimate as follows. Fix e > 0. First we split of the integral J t _ Q which is 
O(e) uniformly in n; so for the rest we have t > e. The last sum of integrals is 
now estimated as follows: Consider first the set of all 9 such that \ct(ip(t, |^ + 

0)> ft + e )\ < £ which 

is a countable disjoint union of open intervals. There 
we get the estimate 0(l).nAn 2 .e 2 (9 3 /3)f ~^ 2n = 0(e), uniformly in n. On the 
complementary set of all 9 where \c t ((p(t, + §f +#)| > £ we use also t > e and 
estimate by 0{l).nAn 2 .^.(9 3 /3)|^ /2n = O(^). The other sums of integrals 
can be estimated similarly, thus L hor (c) goes to for n — > oo. It is clear that 
one can approximate ip by a smooth function without changing the estimates 
essentially. □ 



3.11. Non-smooth curves in the completion of Bi. We have seen in 3.6 
that the completion of Bi in the metric G A lies in the space of Lipschitz maps 
c : S* 1 — > R 2 mod monotone correspondences, that is, rectifiable Frechet immersed 
curves. But how big is it really? We cannot answer this, but we show, in this 
section, that certain non-smooth curves are in the completion. To be precise, if 
c is rectifiable, then we can assume c is parametrized at constant speed \c$\ = 
L/2ir where L is the length of the curve. Therefore c e = (L/2Tr)e ia{ ^ for some 
measurable function a(9) giving the orientation of the tangent line at almost every 
point. We will say that a rectifiable curve c is 1-BVif the function a is of bounded 
variation. Note that this means that the derivative of a exists as a finite signed 
measure, hence the curvature of c - which is (2ir/L)a' - is also a finite signed 
measure. In particular, there are a countable set of 'vertices' on such a curve, 
points where a has a discontinuity and the measure giving its curvature has an 
atomic component. Note that a has left and right limits everywhere and vertices 
can be assigned angles, namely a + (9) — ct-(6). 

Theorem. All 1-BV rectifiable curves are in the completion of Bi with respect to 
the metric G A . 

Proof. This is proven using the following lemma: 

Lemma. Let c(t,6), < t < 1 be an open path of smooth curves c(t) and let 
a(t,9) = arg(ce(6 l , tj). Assume that 

(1) the length of all curves c(t) is bounded by C\, 

(2) \ct\<C 2 , for all(t,6), 

(3) For all t, the total variation in 9 of a(6, t) is bounded by C3 and 

(4) the curvature of c(t) satisfies \k c ^(0, t)\ < C±/t for all 9. 



Then the length of this path is bounded by Ci(\J~C\ + 2V 'AC3C4). 
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To prove the lemma, let s* be arc length on c(t) and estimate the integral: 



/ 

Jc(t) 



(l + AK(c(t))(t,6) 2 )(c u p-} 2 \c e \d6<Cl(C 1+ A / k\ , t) ds t ) 
\ c e\ Jc(t) 



Jc(t) ds t 



<C 2 (C 1 + A^C 3 ). 

Taking the square root of both sides and integrating from to 1, we get the result. 

We apply this lemma to the simplest possible smoothing of a 1-BV rectifiable 
curve c : 

c(i,0) = ^L- / c o (0 - i P )e" p2/2t2 dip = -^L- / cv{v)e- {e -* )2/2t2 d<p^ < t < 1. 
V27rt Jr V27TWR 

Note that t is the standard deviation of the Gaussian, not the variance. We assume 
Co has a constant speed parametrization and c' = (i/27r)e M as above, where a' is 
a finite signed measure. Thus: 

L 



C0 



cee 



iL 



ia(e- v )- v y2r 



d<p 



{difi) 



(27r) 3 / 2 i , 

Moreover, using the second expression for the convolution and the heat equation 
for the Gaussian, we see that q = tcgg. We now estimate: 

| col < L/2tt, hence length(C t ) < L 

< e -*- nL fl 2t2 \a'\{d V ) 



I cee | < 



(2' 

< sup (J2 



e fx-^V2t 2 ^-Var(arg(c )) = 



\d0 < 



L 
27 



2nt 

\ct\=t\c e e\=0(l). 



2 /* 2 dO)(j W{dip)\) = ^Var(arg(c )) 



To finish the proof, all we need to do is get a lower bound on |ce|. However, 
|ce| can be very small if the curve c has corners with small angles. In fact, c can 
even double back on itself, giving a 'corner' with angle 7r. We need to treat this 
as a special case. When all the vertex angles of c are less than n, we can get a 
lower bound for Icel as follows. We start with the estimate: 



\ce{0)\ = 
> 



L- / e^ e -^ 2 / 2t2 dip/e la ^ 

cos(a(8 -ip)- a{e))e-* 2 ' 2t2 dip 



2nt 
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We breakup the integral over R into 3 intervals (-00, 6-8/2], [6-8/2,6+8/2], [6 + 
8/2, +00) for a suitable 8. If t is sufficiently small, the integral of the Gaussian 
over the first and third intervals goes uniformly to and, on the middle interval, 
goes to 1. Thus it suffices to estimate the cos in the middle interval. We use a 
remark on BV functions: 



Lemma. For any BV function f(x) and any C > 0, there is a 8 > such that on 
every interval I of length less than 8, either f\j has a single jump of size > C or 
max(/|j) - min(/| 7 ) < C. 

In fact, let C — e be the size of the largest jump in / less than C and break up 
the domain of / into intervals Ji on each of which the variation of / is less than 
e/2, big jumps being on their boundaries. If 8 is less than the minimum of the 
lengths of the Ji, we get what we want. 

Now let 7r — [3 be the largest vertex angle of the curve cq- Using the last lemma, 
choose a 8 so that on every interval / in the (9-linc of length less than <5, either / 
contains a single vertex with exterior angle > (3/3 or maxa|/ — mina|/ < (3/3. 
Now if there is no vertex in [6 - 8/2, 6 + 8/2], then \a(6 - ip) - a(6)\ < [3/3 on this 
interval and our lower bound is: 

\c e (6)\ > cos(/3/3) - o(t). 

On the other hand , if there is such a vertex, say at 6, then a varies by at most 
f3/3 in [6 — 8/2, 6), jumps by at most ir — (3 at 6 and then varies by at most (3/3 
on (6, 6 + 8/2]. Assume 6 < 6 (the case 6 > 6 is similar). Then: 



cos(a(6-ip)-a(6)) > 
Thus: 



cos(/3/3), if if e {6- 6,6 + 8/2] 

cos(tt -[3 + (3/3) = - cos(2/3/3), if <p G [6 - 8/2, 6-6) 



\cg{6)\ > i(cos(/3/3) - cos(2/3/3)) - o(t). 

hence, if t is sufficiently small, we get a uniform lower bound on \c$\. Since 
\ K c t \ < |c06>|/|c6>| 2 , we get the required upper bound both on \nc t \ and on the 
variation of ac t , i.e. J sl \nc t \ an d all the requirements of the lemma are satisfied. 

If Co has a vertex with angle it, we need to add an extra argument, cq certainly 
has at most a finite number of such vertices and we can construct a new curve 
by drawing a circle of radius t around each of these vertices and letting c^ be 
the curve which follows c until it hits one of these circles and then replaces the 
vertex with a circuit around the circle: see Figure 3. Each of the curves c^ is in 
the completion of Bi by the previous argument and the path formed by the Cq^ 's 
also has finite length, hence Co is in the completion. We omit the details which 
are straightforward. 



3.12. The energy of a path as 'anisotropic area' of its graph in R 3 . 

Consider a path t c(t, •) in the manifold Imm(S' 1 , R 2 ). It projects to a path ttoc 
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Figure 3. Approximating 1 — BV curves with zero angle vertices 
by curves with positive angle vertices. 

in B^S 1 ,^. 2 ) whose energy is 

,6 

E G A (TTOc) = \ G* (c) {T c 7T.C t ,T c TT.C t ) dt 

J a 

= \ [ b G?(ct,cl)dt=l f j (l + An{cY){ci,ci)\c e \dedt 

Ja Ja JS 1 



(1) 



= ± [" [ (i + AK(c) 2 )(^^ice, { -^-ice)\c e \dedt 
Ja Jsi \ \ce\ 2 \ce\ 2 I 

= \\ I (l + An(c) 2 ){c u ic e ) 2 ^-d9dt 

Ja Js 1 \ C 0\ 



If the path c is horizontal, i.e., it satisfies (c t , eg) = 0. Then (c t ,icg) — \c t \.\cg\ and 
we have 

(2) E%>X(c) = \j ( (1 + An{c) 2 )\c t \ 2 \c e \d9dt 1 (c u cg) = 

Ja JS 1 

which is just the usual energy of c. 

Let c(t, 9) = (x(t, 9), y(t, 9)) be still horizontal and consider the graph 

$(t,6) = (t,x(t,6),y(t,6))€m. 3 . 

We also have \x t yg — xgy t \ = |det(c t ,cg)| = |c t |.|cg| and for the vector product 
*f x $« = {x t y e - x e y t , -y e ,xg), so we get 

|* t x $ | 2 - (x t y e - x g y t f + y 2 + x 2 = {x 2 + y 2 ){x 2 + y 2 t + 1) = |c e | 2 (|c t | 2 + 1). 

We express now E hoT (c) as an integral over the immersed surface ScK 3 param- 
eterized by $ in terms of the surface area d/is = \&t x &e \ d9 dt as follows: 

*&{c) = \\ I (l + ^(c) 2 ) ' Ct| ' |Cg| A * t X <Pg\d9dt 
Ja JS 1 \ 9 t X 9g\ 

2 



= \ I (1 + Ak{c) 2 ) Ct d^ s 
JWMxS 1 ^\ct\ 2 + l 
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Next we want to express the integrand as as a function 7 of the unit normal 
ns = ($t x $ e )/|$ t x Let e = (1,0,0), then the absolute value of the 

^-component 4 of the unit normal ns is 



{eo,ns)\= 7= = and 



Thus for horizontal c (i.e., with c t -Lcg) we have 



Horizontal energy as anisotropic area. 



(3) 



E h G °I(c) = \\ (1 + A K (cf)^pM==du 
JK^xs 1 v 1 H4I 



Here the final expression is only in terms of the surface S and does not depend 
on the curve c being horizontal. This anisotropic area functional has to be mini- 
mized in order to prove that geodesies exists between arbitrary curves (of the same 
degree) in B i (S 1 ,R 2 ). Thus we are led to the 

Question. For immersions co,ci : S 1 — > M 2 does there exist an immersed surface 
S = (ins[ .i], c) : [0, 1] x S 1 — > R x R 2 such that the functional (3) is critical at S? 

A first step is: 

Bounding the area. For any path [a, i] 9 f h c(t, ) the area of the graph 
surface S = S(c) is bounded as follows: 

(4) Area(S)= [ d^s < 2E%°X(c) + max£(c(t, ))(b - a) 

J[a,fc]xS! * 



Proof. Writing the unit normal n s — (4' n S' n l) G S* 2 according to the coordi- 
nates (t, x, y) we have 

I4i + 141 + |n ° s|2 > 141 2 + 14 1 2 + 141 2 = 1 

V 1 l n sl 

Since \n s \d/j,s is the area element of the projection of S onto the (t, y)-plane we 
have 

Area(5)=/ d^ s < [ (1 + Ak(c) 2 ) fe| + |4| + j n ° sl \ W 
< 2£7^(c) + max e(c(t, ))(& - a). □ 
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4. Geodesic equations and sectional curvatures 

4.1. Geodesies on Imm(5 1 , R 2 ). The energy of a curve t \— > c(t, •) in the space 
Imm(5 1 ,R 2 ) is 

E G A(c) = h[ I (l + A K 2 c )(c u c t )\c e \dedt. 
By calculating its first variation, we get the equation for a geodesic: 
Geodesic Equation. 

(1) 



((l + ^) M . Ct ) r (^.^. Ce + A 



-I + Ak 2 \c t \ 2 , .Wctl 2 ), . 

• ice 



Proof. From 2.2 we have 



and 



, A _ (ic sg ,cge) (ic e ,c s gg) (c s g,cg) 
K(C)s ~ M 3 + M 3 M 2 ■ 



„„,C0) (cge,ic g ) . 

M \ c e\ 
\ c e\e . / \ I I . 

-C0 + K(c)|Ce|?C e . 



Now we compute 



d s \ E(c) = ia s | / / (l + AK 2 )(c t ,c t )M^ctt 
= T f (A« Ks | Ce || Ct | 2 + (l + ^ 2 )( Cst , c 0| Ce | + ii^|c t | 2 %^) 



\cA 2 \c t \ 2 \c t \ 2 

AK(ic se ,Cgg}- -k + An{ice, c s ee) -. — - 3Ak 2 {c s 0, Cg)- — - 

la Js^ K \ce\ M 2 \ce\ 

[c., ((1 + An 2 )\c B \c t ) t + (l±^!M! C0 ) e )) dm 
c„ ((1 + An 2 )\c 6 \c t ) t + (i±^!^! ce ) )) Mdt 



J ^c s ,-((l + A K 2 )\ C e\c t ) t + F )dedt 



where 



jp a l c *l 2 • ■ a i i i2n ic e , i , 2 Me* c | c* | 2 . 
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2 1 + An 2 \c t \ 2 



-ce 



\c s \ 2 \c e \ 

Substituting the expression for cee and simplifying, this reduces to 



F = 



-1 + Ak 2 \c t \ 



c e + A( K \c t \ 2 )e 



ic e 
\ce\ 2 



2 \cg 

which gives the required formula for geodesies. 

Putting A = in 4.1.1 we get the geodesic equation for the _ff°-metric on 
Imm(S\R 2 ) 



(2) 



(l^)« = -2(n^r) 



4.2 Geodesies on B i (S 1 , R 2 ). We may also restrict to geodesies which are per- 
pendicular to the orbits of Diff(S' 1 ), i.e. (ct,cg) = 0, obtaining the geodesies in 
the quotient space B^S 1 , R 2 ). To write this in the simplest way, we introduce the 
'velocity' a by setting c t = iace/\ce\ (so that \c t \ 2 = a 2 ). When we substitute this 
into the above geodesic equation, the equation splits into a multiple of ce and a 
multiple of ice- The former vanishes identically and the latter gives: 



{(l + An 2 )\c e \a) 



ice 



Ak< 



{{l + An 2 )\c e \a) t 



2 

-1 + Ak 2 



\('o 



CO 



ice 

M 



or 



n\ce\a 2 + A( 



(na 2 )b 
\ce\ 



If we use derivatives with respect to arclength instead of and write these with 
the subscript s, so that f s = ,fe/\ce\, this simplifies. We need: 



\cg\t 



(ce,c t e) _ (ce9,c t ) 



{cee, ice) 



-an\ce\ 



\ce\ \ce\ |car 

which gives us a simple form for the equation for geodesies on _B i (S' 1 ,R 2 ): 



(1) 



((l + A K 2 )a) t = 



1 -\- 3/4 Ax 
t ^ Ka + A ( Ka > s 



Finally, we may expand the ^-derivatives on the left hand side, using the formula 
K t = a,K 2 +a ss noted in 2.2.7; we also collect all constraint equations that we chose 
along the way: 



(2) 







a t = 



(c t ,c s ), c t = aic Sl k = (c ss ,ic s ) 
\k,cl 2 + A(cl 2 (k, ss - i« 3 ) + An s aa s + 2«a 2 ) 
1 + Ak 2 



Handle this with care: Going to unit speed parametrization (so that f s is really 
a holonomic partial derivative) destroys the first constraint 'horizontality'. This 
should be seen as a gauge fixing. 
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12 2 1 f a d\ 

(1) St — —ans, at — -kna , Kt = an H — — 

z s \ s J I 



4.3. Geodesies on Bi(S 1 , R 2 ) for ^4 = 0. Let us now set A = 0. We keep looking 
at horizontal geodesies, so that (ct,ce) — and c t = iac6»/|ce| for a e C 00 (S' 1 ). We 
use the functions a, s = |cg|, and n. We use equations from 4.2 but we do not use 
the anholonomic derivative: 

2 aeg a^se 
S \ S I 8 s z s^ 

We may assume that s\ t =o is constant. Let v(9) = a(0,8) be the initial value for 
a. Then from equations (1) we get 

= -an = -2^ =► log( S a 2 ) 4 = 
s a 

so that sa 2 is constant in t, 

(2) s(i, 6>)a(t, 6>) 2 = s(0, <9)a(0, 6>) 2 = v(6>) 2 , 

a smooth family of conserved quantities along the geodesic. This leads to the 
substitutions 

S — — ~, K — Z— ^ 

a z a z 
which transform the last equation (1) to 

a\ A a t a&a ee , a 6 a v g a b a 2 e 

(3) att-4 tt^-^ 5 4- = °> a{0,6) = v{6), 



a nonlinear hyperbolic second order equation. Note that (2) implies that wherever 
v = then also a = for all t. For that reason, let us transform equation (3) into 
a less singular form by substituting a = vb. Note that b = 1/ ^Js. The outcome is 

(4) (b-% t = - V 4(b 3 ) ee 2vv e (b 3 )e ^6 3 , b(0,8) = 1. 



4.4. The induced metric on Bij(S ,M ) in a chart. We also want to compute 
the curvature of B^S 1 ,^ 2 ) in this metric. For this, we need second derivatives 
and the most convenient way to calculate these seems to be to use a local chart. 
Consider the smooth principal bundle ir : Imm/(S' 1 , R 2 ) —> Bij(S 1 ,M. 2 ) with 
structure group Diff(S' 1 ) described in 2.4.3. We shall describe the metric in the 
following chart near C € Bjj(S' 1 ,R 2 ): Let c € Immj(S' 1 ,M 2 ) be parametrized by 
arclength with ir(c) = C of length L, with unit normal n c . We assume that the 
parameter 8 runs in the scaled circle below. As in the proof of 2.4.3 we consider 
the mapping 

V : C°°(Sl (-e, e)) - Imm/OSi, M 2 ), Q(c) := 4>(C°°(Sl (s, e))) 

= c(6) + f(8)n c (8) = c(8) + f{6)id{6), 

Tro4>:C°°(Sl(-e,e))^B hf (S\R 2 ), 

where e is so small that ip(f) is an embedding for each /. By 2.4.3 the mapping 
(it o ip)" 1 is a smooth chart on Bij(S^, R 2 ). Note that: 

VK/)' = c + fie' + fid' = (1 - fn c )c + f'ic 

W)" = c" + f'ic' + 2 fid' + fid" = -(2/'k c + K)c' + (« c + /" - /« 2 ) iC ' 
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1 



- f^ c ? + I' 2 



([l-f Kc )ic'-f'c'), 



Tfip.h = h.ic' e C°°(S' 1 ,M 2 ) = T W) Imm f (S'i,R 2 ) 
HI - fn c ) hf 

V(i-/« c ) 2 + /' 2 (i-/k c ) 2 +/' 



(T^.ft^ = -jML fKc) n W) 

7(1 - /«c) 2 + J' 2 



</-(/)' 



K W) = ((1_ /Kc )2 + ^2)3/2 W{f)'Mf)") 



= « c + r - 2/zcg - ff" Kc + / 2 ^g + 2/'^ c + //X 

((i - fn c y + /' 2 )3/2 

Let G A denote also the induced metric on Bij(S\, R 2 ). Since 7r is a Riemannian 
submersion, T w) -n : (A^ (/) , G^ (/) ) -> {Bij{S\, R 2 ), G^ (/)) ) is an isometry. 
Then we compute for / e C°°(5i, (-e,e)) and ft, fc € C^S^R) 

((tt o WG A ) f {h, fc) = G^ (/)) (^(tt o ^)ft, T f (n o ^)k) 
= G^ f) ((T f i ) .h) ± ,(T f i,.k) ± ) 

= J (1 + A4 (/) )((T / ^.ft)- L , (T / V.fc)- L )|V(/)'l ^ 



) 2 + /' 2 



This is the expression from which we have to compute the geodesic equation in 
the chart on B u (Sl, R 2 ). 

4.5. Computing the Christoffel symbols in B i j(S} j ,R 2 ) at C = ir(c). We 
have to compute second derivatives in / of the expression of the metric in 4.2. For 
that we expand the two main contributing expressions in / to order 2, where we 
put k = K c . 

K 4>(f) = 

= (1 - 2/k + / 2 « 2 + /' 2 )- 3 / 2 ( K + /" - 2/« 2 - f f" k + / 2 « 3 + 2/' 2 « + //V) 
= « + (/" + /^ 2 ) + (/ 2 « 3 + |/' 2 « + //'«' + 2//"k) + 0(/ 3 ) 
(1 - fnf(l - 2fn + /V + f' 2 )- 1 / 2 = 1 /ft if' 2 + 0(/ 3 ) 
Thus 



(1 + At 2 ) (1 /Kc)2 = 1 + Ak 2 + 2A/"« + A/« 3 - /«- 

/ (1-/Kc) 2 + /' 2 
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- if' 2 + Afn 4 + A\fn 2 + 2A//W + Af" 2 + 4A.ff" K 2 

and finally 

(1) Gf(h,k) = ((TTOlP)*G A ) f (h 7 k) = 

= f hk ((1 + An 2 ) + (2Af"n + Afn 3 - /«) + -\f 

+ A(4//V + /V + ±.f V + 2//W + f" 2 ) + C>(/ 3 )) d6. 
We differentiate the metric 

dG A (f)(l)(h,k) = [ hk(2Al"n+ (An 3 - n)l + 4Alf"n 2 +4Afl"n 2 + 
J si v 

+ 2Afln 4 + (An 2 - l)f'l' + 2Alf' kk! + 2AH'nn! + 2Af"l" + 0(/ 2 )) d9 
and compute the Christoffel symbol 

- 2Gf(T f (h, k), I) = -dG A (.f)(l)(h, k) + dG A (f)(h)(k, I) + dG A (.f)(k)(l, h) 

= [ l((An 3 - K + 2Ann'f + AAn 2 f" + 2An i f)kh 
J si v 

+ (2 An + AAk 2 f + 2Af")(h"k + hk") 

+ [An 2 ]' - f + 2A K K'f)(h'k + hk') + 0(/ 2 )) dff 

- [ (l'{An 2 f'hk- f'hk + 2Ann' fhk) 

+ l"(2Anhk + 4An 2 fhk + 2Af"hk) + 0(/ 2 )) d9 

= [ l((An 3 - K-2AK")hk-4AK'{h'k + hk')-4AKh'k' 
J si v 

+ (-2Af^-f" + 2Ak 4 / - 6Ak' 2 / - 6Akk" f - WAndf + An 2 f")hk 

- (2f + AAf" + \2Akk'! + 6An 2 f')(h'k + hk') 

- 2(4A K 2 f + 2Af")h'k' + 0(/ 2 )) d6. 

Thus 

Gf(T f (h,k),l) = 

= ( i({\k- \Ak 3 + AK")hk + 2An'{h'k + hk') + 2AKh'k' 
J si v 

- (-A/ (4) - if" + An A f - 3An' 2 f - 3Akk" f - 5Akk' f + \Ak 2 f")hk 
+ (/' + 2Af" + QAnn'f + 3A K 2 f')(h'k + hk') 

+ (AAk 2 / + 2Af")h'k' + 0(/ 2 )) d0. 
At the center of the chart, for / = 0, we get 
G A (T (h,k),l) = 
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= / - \ Ak5 + An")hk + 2An'{h'k + hk 1 ) + 2Anh'k'^ 

■ {\k - \Ak> + An")hk + 2An'{tik + hk') + 2Anh'k' 



1 + An 



2 



yi + An 2 ) d9 



= G, 
so that 



a ((\k- \-Ak 3 + An")hk + 2An'(h!k + hk') + 2Anh'k' 
\ l + An 2 ' 



(2) T (h,k) 



(i K _ ^Ak 3 + An")hk + 2An'(h'k + hk 1 ) + 2Anh'k' 



1 + An 2 

Letting h = k = f t , this leads to the geodesic equation, valid at / = 0: 
_ ( i K _ t;Ak 3 + An")f 2 + W/t/ t ' + ^Mft? 
Ht 1 + Ak 2 

If we substitute a for f t and a t for / tt , this is the same as the previous geodesic 
equation derived in 4.2 by variational methods. There is a subtle point here, 
however: why is it ok to identify the second derivatives at and ftt with each other? 
To check this let c{9) + (£ai(0) + ^a 2 (6))ic'(0) be a 2-jet in our chart. Then if we 
reparametrize the nearby curves by substituting 9 — ^-aia[ for 9, letting 

2 2 t 2 2 2 

c(£, 9) = c(9 - ^a^) + (iai(0 - ^a^a'J + —a 2 {9 - £aiai))ic(0 - ^Wi)' 

2 t 2 

= c(9) - ( t T a 1 a' 1 )c'{9) + (iai(0) + -a 2 {9))ic' (9) mod £ 3 

then (c',c t ) = mod £ , hence this 2-jet is horizontal and {ctt,ic') = 02 mod £ 
as required. 

4.6. Computation of the sectional curvature in E>ij{S\, M. 2 ) at C. We now 
go further. We use the following formula which is valid in a chart: 

(1) 2R f (m,h,m,h) = 2Gf(R f (m,h)m,h) = 

= -2d 2 G A (f)(m,h)(h,m) + d 2 G A (f)(m,m)(h,h) + d 2 G A (f)(h, h)(m,m) 

- 2G A (T(h, to), r(m, h)) + 2G A (r(m, m),T(h, h)) 

The sectional curvature at the two-dimensional subspace Pf(m, h) of the tangent 
space which is spanned by m and h is then given by: 

Gf(R(m,h)m,h) 

(2) k f (P(m,h)) = - Mmi2 _ G A {m , hr 

We compute this directly for / = 0. From the expansion up to order 2 of G A (h, k) 
in 4.5.1 we get: 

(3) ^d 2 G A {0){m,l)(h,k)= j hk(-\m'l'+ 



si 



+ A(2(ml" + m"l)n 2 + mln A + \m'l'n 2 + {ml 1 + m'l)™' + m"l"^j 



d9 
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Thus we have: 

-d 2 G A (0)(m,h)(h,m) + ±d 2 G A (0)(m,m)(h,h) + ±d 2 G A (0)(h, h)(m,m) = 
= —2 J ^ hm (—^m'h l + 

+ A(2(mh" + m"h)n 2 + mhn 4 + \m!h!n 2 + (mh' + m'h) kk' + m"/i")) d6 
+ [ hh (-\m' 2 + A[^mm!'n 2 + m 2 k 4 + \m! 2 k 2 + 2mm! kk! + m" 2 )) dB 

+ J mm (-\titi + A(4hh"n 2 + hhn A + \titin 2 + 2hti kk' + h"h"^ 

= J (j{Ak 2 - l){mti - m'h) 2 + A(mh" - m"/i) 2 ) d6. 

For the second part of the curvature we have 

- G {T {h, m), r (m, h)) + G (T (m, m), T (h, h)) = 



d6 



■■ J ^ \ Ak z + AK")hm + 2AK'(h'm + m'h) + 2AKh'm') - 

+ / {{\k - \Ak 3 + AK")m 2 + AAK'mm' + 2Akto' 2 ) 
■' s l 

_ i Ak 3 + AK")h 2 + AAK'hti + 2AKh' 2 ^ - d0 

- [ ((Ak 2 -A 2 k 4 + 2A 2 kk" -AA 2 K l2 ){mh' -m'h) 2 )- 
Jsl K > 1 



+ Ak 2 



Ak 2 
d9 



+ Ak 2 
Thus we get 

Ro{m, h, m, h) = G A (Ro(m, h)m, h) = 

= J (\{Ak 2 - l)(mh' - m'h) 2 + A{mh" - m"hf) d9 

+ J i ({Ak 2 - A 2 k a + 2A 2 kk" - AA 2 K l2 )(mh' - m'h) 2 ) 
Letting W = mh' — hm' be the Wronskian of m and h and simplifying, we have: 



(4) 



Ro(m, h, m, h) — 

„2 i\2 , a A2 // o A2,j2 s 



Js! 2(1 + A,?) I J s! 



What does this formula say? First of all, if supp(m) n supp(/i) = 0, the sectional 
curvature in the plane spanned by m and h is 0. Secondly, we can divide the curve 
c into two parts: 



c\ = set of points where kk" < 2(k') 2 + 

= set of points where kk" > 2{k') 2 + ( A ~ 1 ^~ K 
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Note that if A is sufficiently small, = and even if A is large, c~ A need not 
be non-empty. But if supp(m), supp(/i) C c A , the sectional curvature is always 
negative. The interesting case is when supp(m), supp(ft-) C c\. We may introduce 
the self-adjoint differential operator on L 2 (S' 1 ): 

(An 2 - If - AA 2 nn" + 8AV 2 
1 1 2A{l + An 2 ) ■' 

so that R = — A(SW, W). The eigenvalues of S tend to — oo, hence S has a 
finite number of positive eigenvalues. If we take, for example, m = 1 and h 
such that b! is in the span of the positive eigenvalues, the corresponding sectional 
curvature will be positive. In general, the condition that the sectional curvature 
be positive is that the Wronskian W have a sufficiently large component in the 
positive eigenspace of S. The special case where c is the unit circle may clarify 
the picture: then 

Sf = f " + 2A(1 + A) f 
and the eigenfunctions are linear combinations of sine's and cosine's. It is easy 
to see that for any A, a plane spanned by m and h of pure frequencies k and I 
will have positive curvature if and only if k and I are sufficiently near each other 
(asymptotically \k — l\ < \A— l\/\/A + a 2 ), hence 'beat' at a low frequency. 

4.7. The sectional curvature for the induced i/°-metric on Bij^S^M. 2 ) 
in a chart. In the setting of 4.2 we have for / e C°°{S\ j , (— e,e)) and h, k <G 
C°°(Sl,R) 

(1) G° f (h, k) = ((tt o ipyG°) f (h, k) = G% (f)) (T f (n o IXtt o ^)k 

= G° (/) ((T / ^) ± ,(T / M) ± ) 
hk(l - fn c ) 2 



: (19 



At the center of the chart described in 4.4, i.e., for / = 0, the Christoffel symbol 
4.5.2 for A = becomes 

(2) T (h,k) = \n c hk 
The curvature 4.6.4 at / = for A = becomes 

Ro(m, h, m, h) = Go(Ro(m, h)m, h) = 

(3) =-\\ {h'm-hm') 2 d6 = -\ [ W{m,h) 2 d9 

J sl Jsi 

and the sectional curvature ko(P(m,h)) from 4.5.2 for A = and / = is non- 
negative. 



In the full chart 4.2, starting from the metric 4.6.1, we managed to compute 
the full geodesic equation not just for / = but for general /, so long as A = 0. 
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The outcome is 



Kc h 2 -\k c {1- fK c )h 2 + {\h?f" + 2hh'f) 



(4) 



1 - /«c ((1 - fn c ) 2 + f' 2 ) 

Kc h 2 r 2 , i Kc (i-fK C )h 2 r 2 -ih 2 f 2 f" 



2 



(l-/«c)((l-/«c) 2 + /' 2 ) ((l-/«c) 2 + /' 2 ) 

The geodesic equation is thus 

1 - /«c ((1 - fKc? + / 6 2 ) 

K J 2 f t h 2 f « c (l - /K c )/ t V - |/ 2 / e 2 J 



C5) j 

(l-/«c)((l-/«c) 2 + /fl 2 ) ((l-/«c) 2 + /fl 2 ) 2 

For A > we were unable to get the analogous result. 

5. Examples and numerical results 

5.1. The geodesies running through concentric circles. The simplest pos- 
sible geodesic in Bi is given by the set of all circles with common center. Let C r be 
the circle of radius r with center the origin. Consider the path of such circles C r / t ) 
given by the parametrization c(t,6) = r(t)e l6 , where r(t) is a smooth increasing 
function r : [0, 1] — > K>o- Then n c (t, 9) = l/r{t). If we vary r then the horizontal 
energy and the variation of this curve are 

E^(c) = \j o j s (l + A/r 2 y t rd6dt 

dsU^Kc) = I' I (l + *)r s (-r tt - £zA^) rdBdt 
so that c is a geodesic if and only if 

Also the geodesic equation 4.1.1 reduces to (1) for c of this form. 



The solution of (1) can be written in terms of the inverse of a complete elliptic 
integral of the second kind. More important is to look at what happens for small 
and large r. Asr^ 0, the ODE reduces to: 

r 2 
2r 

whose general solution is r(t) = C(t — t ) 2 for some contants C, t . In other words, 
at one end, the path ends in finite time with the circles imploding at their common 
center. Note that r' — > as r — > but not fast enough to prevent the collapse. 
On the other hand, as r — > oo, the ODE becomes: 
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whose general solution is r(t) = C(t — to) 2 / 3 for some constants C, to- Thus at the 
other end of the geodesic, the circles expand forever but with decreasing speed. 

An interesting point is that this geodesic has conjugate points on it, so that it 
is a extremal path but not a local minimum for length over all intervals. This is a 
concrete reflection of the collapse of the metric when A = 0. To work this out, take 
any f(9) such that ^ fd9 — and any function a(t). Then X = f(8)a(t)d/dr 
is a vector field along the geodesic, i.e. a family of tangent vectors to B e at each 
circle C r (t) normal to the tangent to the geodesic. Its length is easily seen to be: 

\\X\\ 2 c rW = {r{t) + ^)a{tf f{6) 2 d6. 

We need to work out its covariant derivative: 

V* t (X) = m a t l+r Cr (r t l, m al). 

Using a formula for the Christoffel symbol which we get from 4.2.2 by polarizing, 
and noting that n = 1/r, n s = 0, we get: 

d .... / 1 -A/r 2 \ d 



V* t (X) = f(6)a t -+f(e)ar t Q 



(r + A/r)J dr 



f(9)(r + A/r) 1/2 ((r + A/r) 1/2 a) 



1/2 J\ d 
dr 



(This formula also follows from the fact that the vectors (r + A/r)~ x l 2 d /dr have 
length independent of t, hence covariant derivative zero.) Jacobi's equation is 
therefore: 

(2) m (r + A/r) ^ ((r + A/r) 1/2 a) J- + R(X, r t ^) (r t ±) = 0, 

where R is the curvature tensor. For later purposes, it is convenient to write this 
eqation using r as the independent variable along the geodesic rather than t and 
think of a as a function of r. Note that for any function b along the geodesic, 
b t = b r r t and 

(1 — A/r 2 ) 

ht = b rr r 2 + b r r tt = (b rr ~ —, — , , / A ) r t ■ 

2{r + A/r) 

Then a somewhat lengthy bit of algebra shows that: 

1 / Is 1 / 1 



(r + A/r) 2 ((r + A/r) 2 a) = (r + A/r) 4 ( (r + A/r) 4 a) r 2 + F(r)ar 

\ J tt V / rr 



5 /1-^A 2 \ 2 
where F(r) = - — [ — — {-j-} + 



16 V r + A/r) 2r 3 (r + A/r)' 

To work out the structure of R in this case, use the fact that the circles C r 
and the vector field d/dr are invariant under rotations. This means that the map 
/ i ► R(d/dr, f d/dr) (d/dr) has the two properties: it commutes with rotations 
and it is symmetric. The only such maps are diagonal in the Fourier basis, i.e. 
there are real constants A„ such that 

r>(xix ( cos ( n ^)9/dr\, a/a . (cos(n9)d/dr 
R { d ^ UnUwdrF^ = Hsin(n6)d/dr 
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To evaluate A„, we take the inner product with cos(n#) (or sin(n#) and use our 
calculation of Ro(m, h, m, h) in section 4.6 to show: 

R (<j- , CO s(n0) ^) (^), cos(ne) ) = R , cos(n#) ^ , ^ , cos(nfl) ^ 



where 



... , (i , sm(n0) . d Jlr 2 cos(n6>) 
W = 1 .— cos(nfl) = -n — - — - and W = — W = -n 2 V^- 

as r as r 2 



Simplifying, this gives: 

A„|| cosM)|-|| 2 = J^- ^-f^ n^iuO) + ^cos 2 (n0)) 



(l-A/r 2 ) 2 2 , A 



, , n 2 ir + — n 4 7r 
2(r + A/r) r 3 

hence 

(1 - A/r 2 ) 2 2 A 4 

2(r + ^L/r) 2 r 3 (r + A/r) ' 

Thus for X = cos(n9)a n (t)d/dr, if we combine everything, Jacobi's equation reads: 

(r + A/r)~* ((r + A/r)*a n ^ 



(3) 



/_ (l-A/r 2 ) 2 A _ n 

I 2(r + A/r) 21 sJ + rZ(r + A/r) [ 2> ) 



Calling the right hand side the potential of Jacobi's equation, we can check that 
for each n, the potential is positive for small r, negative for large r and it has one 
zero, approximately at y/2An for large n. Thus, for small r, these perturbations 
diverge from the geodesic of circles. For large r, if we write b n = (r + Ajr) X ^a n , 
then Jacobi's equation approaches: 

_ n 2 - 0.625 ^ 

\O n ) rr ~ ^ ^ °n- 

This is solved by b n — cx x + c'x x where A, A' are solutions of A 2 — A = — (n 2 — 
0.625)/2. For n = 1, A, A' are real and 6„ has no zeros; but for n > 1, A, A' have 
an imaginary part, say i^ n , and 

b n « v / 7"(ccos(7„ log(r)) + c' sin(7 log(r))) 

with infinitely many zeros. 

Figure 4 shows the solution for n = 3 which approaches as r — > 0. The first 
zero of this solution is about 10.77\/A, making it a conjugate point of r = 0. For 
other n, the first such conjugate point appears to be bigger, so we conclude: on any 
segment < n < r 2 < 10.77\/A the geodesic of circles is locally (and presumably 
globally) minimizing. 
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Figure 4. The potential in the Jacobi ODE and its solution for 
an infinitesimal triangular perturbation of the circles in the geo- 
desic of concentric circles. Note the first conjugate point at 
10.77%/A 

5.2. The geodesic connecting two distant curves. For any two distant curves 
C\, C2, one can construct paths from one to the other by (a) changing G\ to some 
auxiliary curve D near C\ , (b) translating D without modifying it to a point near 
C2 and (c) changing the translated curve D to C2. If C\ and C2 are very far 
from each other, the energy of the translation will dominate the energy required 
to modify them both to D. Thus we expect that a geodesic between distant curves 
will aymptotically utilize a curve D which is optimized for least energy translation. 
To find such curves D, heuristically we may argue that it should be a curve such 
that the path given by all its translates in a fixed direction is a geodesic. 

Such geodesies can be found as special cases of the general geodesic. We fix 
e = (1,0) as the direction of translation and assume that the path {D + te] is a 
geodesic. We need to express this geodesic up to order 0{t 2 ) in the chart used in 
section 4.4. Let c(s) be arc length parametrization of D and 6(s) be the orientation 
of D at point c(s), i.e. c s = cos(6*) + ism(9). Then a little calculation shows that 
if we reparametrize nearby curves via s = s — (e, c s )t, then the path of translates 
in direction e is just: 
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Thus, in the notation of 4.2, a — — sin(0), hence a s = — cos(9)k and, moreover, 
at = cos 2 (8)k. Substituting this in the geodesic formula 4.2.1, we get 

(1 + An 2 )cos 2 (9)n = 

= 1^M+A({k ss - y)sin 2 (6») +4cos(6»)sin(6l) KKs +2k 3 cos 2 (6I)). 

Since n = 9 S , this becomes, after some manipulation, a singular third order equa- 
tion for 0(s): 

9 SSS = Acot(9)9 s 9 ss + (I - cot 2 (9))9 S (9 2 S - \). 

One solution of this equation is 8(s) = -^=, i.e. a circle of radius VA. In fact, this 
seems to be the only simple closed curve which solves this equation. However, if 
we drop smoothness, a weak solution of this equation is given by the C , piecewise 
C 2 -curve made up of 2 semi-circles of radius \I~A~ joined by 2 straight line segments 
parallel to the vector e and separated by the distance 2\J~A (as in figure 5). Note 
that such 'cigar'-shaped curves can be made with line segments of any length. 

1.5 1 i i i ii 




Figure 5. On the top, the geodesic joining circles of radius 1 at 
distance 3 apart with A = . 1 (using 20 time samples and a 40-gon 
for the circle). On the bottom, the geodesic joining 2 'random' 
shapes of size about 1 at distance 5 apart with A — .25 (using 20 
time samples and a 48-gon approximation for all curves). In both 
cases the middle curve which is highlighted. 

A numerical approach to minimize E^{c) for variations c with initial and end 
curves circles at a certain distance produced the two such geodesies shown in Figure 
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Figure 6. The forward integration of the geodesic equation when 
A = 0, starting from a straight line in the direction given by a 
smooth bump-like vector field. Note that two corner like singu- 
larities with curvature going to oo are about to form. 

5. Note that the middle curve is indeed close to such a 'cigar'-shape. However, the 
width of this shape is somewhat greater than 2\J~A: this is presumably because 
the endcurves of this path arc not sufficiently far apart. Thus experiments as well 
as the theory suggest strongly that geodesies joining any two curves sufficiently 
far apart compared to their size asymptotically approach a constant 'cigar'-shaped 
(^-intermediate curve made up of 2 semi-circles of radius \J~A~ and 2 parallel line 
segments. We conjecture that this is true. 

5.3 The growth of a 'bump' on a straight line, when A = 0. We have seen 
above that the geodesic spray is locally well-defined when A = 0. To understand 
this spray and see whether it appears to have global solutions, we take that the 
initial curve contains a segment with curvature identically zero, i.e. contains a line 
segment, and that the initial velocity a is set to a smooth function with compact 
support contained in this segment. For simplicity, we take the velocity a to be a 
cubic B-spline, i.e. a piecewise cubic which is C 2 with 5 non-C 3 knots approximat- 
ing a Gaussian blip. The result of integrating is shown in Figure 6. Note several 
things: first, where the curvature is zero, the curve moves with constant veloc- 
ity if we follow the orthogonal trajectories. Secondly, where the curve is moving 
opposite to its curvature (like an expanding circle, the part in the middle), it is 
deccelerating; but where it is moving with its curvature (like a contracting circle, 
the parts on the 2 ends), it is accelerating. This acceleration in the 2 ends, creates 
higher and higher curvature until a corner forms. In the figure, the simulation is 
stopped just before the curvature explodes. In the middle, the curve appears to 
be getting more and more circular. As the corners form, the curve is approaching 
the boundary of our space. Perhaps, with the right entropy condition, one can 
prolong the solution past the corners with a suitable piecewise (7 -curve. 

Although this calculation assumes A = 0, one will find very similar geodesies 
when A is much smaller than l/n 2 1 1/(k s log(a) s ) and k/k ss , so that the dominant 
terms in the geodesic equation are those without an A. In other words, geodesies 
between large smooth curves are basically the same as those with A = 0. 
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5.4 Several geodesic triangles in B e . We have examined dilations, translations 
and the evolution of blips. We look next at rotations. To get a pure rotational 
situation, we consider ellipses centered at (0,0) with the same eccentricity 3 and 
maximum radius 1, but differently oriented. We take 3 such ellipses, with orien- 
tations differing by 60° and 120° degrees. Joining each pair by a geodesic, we get 
a triangle in B e . 

We wanted to examine whether along the geodesic joining 2 such ellipses (a) one 
ellipse rotates into the other or (b) the initial ellipse shrinks towards a circle, while 
the final ellipse grows, independently of one another, ft turns out that, depending 
on the value of A, both can happen. Note that we get similar geodesies by either 
changing A or making the ellipses smaller or larger with A held fixed. For each 
A, we get an absolute distance scale with unit l/VA and, if the ellipses are bigger 
than this, (b) dominates, while, if smaller, (a) dominates. 

The results are shown in Figure 7. We have taken the three values A = 1,0.1 
and 0.01. For each value, on the top, we show the geodesic joining 2 of the ellipses 
as a sequence of curves in their common ambient R 2 . Below this, we show the triple 
of geodesies as a triangle, by displaying the intermediate curves as small shapes 
along lines joining the ellipses. This Euclidean triangle is being used purely for 
display, to indicate that the computed structure is a triangle in B e . Note that for 
A = 1, the intermediate shapes are very close to ellipses, whose axes are rotating; 
while for A = 0.01, the bulges in one ellipse shrink while those of the other grow. 

We can also compute the angles in B e between the sides of this triangle. They 
work out to be 34° when A = 1, i.e. the angle sum for the triangle is 102°, much less 
than 7r radians, showing strong negative sectional curvature in the plane containing 
this triangle. But if A = 0.1 or 0.01, the angle is 77° and 69° respectively, giving 
more than 7r radians in the triangle. Thus the sectional curvature is positive for 
such small values of A. 



5.5 Notes on the numerical simulations. All simulations in this paper were 
carried out in MatLab. The forward integration for the geodesic equation for 
^4 = was carried out by the simplest possible finite difference scheme. This seems 
very stable and reliable. Solving for the geodesies was done using the MatLab 
minimization routine f minunc using both its medium and large scale modes. This, 
however, was quite unstable due to discretization artifacts. A general path between 
two curves was represented by a matrix of points in M 2 , approximating each curve 
by a polygon and sampling the path discretely. The difficulty is that when the 
polygons have very acute angles, the discretization tends to be highly inaccurate 
because of the high curvature localized at one vertex. Initially, in order to minimize 
the number of variables in the problem, we tried to use small numbers of samples 
and higher order accurate discrete approximations to the derivatives. In all these 
attempts, the discrete approximation "cheated" by finding minima to the energy 
of the path with polygons with very small angles. The only way we got around this 
was to use first order accurate expressions for the derivatives and relatively large 
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FIGURE 7. Top Row: Geodesies in three metrics joining the same 
two ellipses. The ellipses have eccentricity 3, the same center and 
are at 60° degree angles to each other. At left, A = 1; in middle 
A = 0.1; on right A = 0.01. Bottom Row: Geodesic triangles in 
B e formed by joining three ellipses at angles 0, 60 and 120 degrees, 
for the same three values of A. Here the intermediate shapes are 
just rotated versions of the geodesic in the top row but are laid 
out on a plane triangle for visualization purposes. 



numbers of samples (e.g. 48 points on each curve, 20 samples along the geodesic, 
hence 2 x 20 x 48 = 1920 variables in the expression for the energy. 

Another problem is that the energy only depends on the path of unparametrized 
curves and is independent of the parametrization. To solve this, we added a term 
to the energy which is minimized by constant speed paramctrizations. This still 
leaves a possibly wandering basepoint, and we added e times another term which 
asked that all points on each curve should move as normally as possible. In 
practice, if the initialization was reasonable, this term was not needed. The final 
discrete energy that was minimized was this. Let Xij be the i th sample point on 
the j th curve Cj. For each estimate the sum of the squared curvature of Cj 

plus the squared acceleration of the parametrization by: 

(The harmonic mean of the segment lengths is used here to further force the 
parametrization to be uniform.) Then, for each (i, j), the four triangles t = {a = 
(i,j),b = (i ± 1, j),c — ± 1)} around are considered and the energy is 
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taken to be: 



n 



((x a - Xb), (x, 



a 



x c )^) 2 + e{(x. 



Xb), (X a 



)(l + _4A(a)). 



i,j,t 



\\x a - X b \\ 



We make no guarantees about the accuracy of this simulation! The results, how- 
ever, seem to be stable and reasonable. 
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